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ABSTRACT 

The  analysis  of  scattered  electromagnetic  radiation  is  examined  in 
terms  of  Mie  theory  establishing  the  detailed  form  for  the  scattering 
functions.  These  parameters  are  used  to  define  the  criterion  whereby 
the  composition  and  number  densities  of  an  atmospheric  aerosol  population 
may  be  assessed.  Laser  radar,  or  LIDAR,  is  the  proposed  remote  sensing 
device  and  the  governing  system  equations  are  developed.  The  problem  of 
data  inversion  is  surveyed  with  emphasis  on  smco tiling  methods,  statistical 
analyses,  and  iterative  techniques.  A  discussion  of  the  numerical 
stability  of  the  solution  is  also  presented.  On  the  basis  of  the  Mie 
model,  inversion  of  the  data  provided  by  the  LI DA?,  probe  is  given  as  a 
rationale  to  effect  an  estimate  of  the  true  particle  number  distribution 
function  in  an  atmospheric  cell. 
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I.  riTRCDL'GTIQi: 


The  advent  of  sophisticated  optical  devices  in  both  the  scientific 


and  defense  communities  has  pointedly  demonstrated  the  need  to  better 


understand  the  primary  propagating  and  interacting  medium  characteristics 


of  the  earth1 s  atmosphere.  The  intent  of  this  paper  is  not  t-  deal  with 


the  purpose  of  such  devices,  -whether  they  are  atmospheric  pollutant  mon¬ 


itors  or  infrared  threat  detectors.  Rather  the  concern  here  is  the  use 


of  these  devices  in  the  real  -world  environment  and  the  limitations  im¬ 


posed  on  them  as  a  result  of  atmospheric  optical  signal  degradation. 


This  degradation  is  the  result  of  manifold  scattering  effects  such  as 


Brillouin,  Rayleigh,  Mie,  Fluorescent  and  Raman  scattering  and  includes 


the  atmospheric  constituents1  differential  absorption  of  optical 


radiation.  00 


The  atmospheric  scattering  and  absorption  centers  include  the  air 


molecules,  suspended  particulates  (aerosols)  and  water  droplets  in  fog, 


rain  or  hail.  The  presence  of  aerosols  strongly  determines  the  trans¬ 


mission  characteristics  of  the  atmosphere  and  may  even  dominate  the 


propagation  characteristics.  See  Table  I  and  figure  1  for  a  comparison 


of  particle  sizes  and  the  range  of  number  densities. 


At  the  small  end  of  the  size  soectrum  are  the  air  molecules  them¬ 


selves  with  radii  of  10~^*  m  and  a  concentration  of  10^  per  cm^.  The 


2  L 

large  end  includes  rain  drops  with  radii  of  10  to  10^*m  and  a  cor.cen- 
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tration  of  10  to  10  per  cm  .  The  "small'’  particles,  defined  by  a 


radius  much  less  than  the  wavelength  of  radiation  used  as  the  probe. 
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Humbar  density  of  particles  dlj/dlof'  r  (cm 


-  jey  the  Rayleigh  scattering  model,  ex] _ = _ _ 

decreases  proportionally  as  the  particle  radius  to  wavelength  ratio  to 
the  fourth  power.  It  is  this  fourth  power  law  for  nolecular  scattering 
that  Rayleigh  used  to  explain  the  blue  color  of  the  daytiae  sky. 

TABLE  I 

-3\ 


Particle  tyoes 

Particle  Radius  Ud 

Concentration 

Aitken  nuclei  5 

10“3  to  10" 2 

102  tc  10' 

haze  particles 

10“ 2  to  1 

10  to  103 

fog  droplets 

1  to  10 

10  to  102 

cloud  droplets 

1  to  10 

10  to  300 

*  Initial  nuclei  present  in  the  atmosphere  iron  the  earth*  s  surface 
fit I  and/or  neteoric  sources. 


For  wavelengths. larger  than  those  of  visible  light,  this  scattering  is 
very  snail.  As  the  particle  sine  approaches  that  of  the  radiation  wave¬ 
length,  the  complete  1-lie  theory  is  required  to  explain  the  scattering 
characteristics.  The  dependence  on  wavelength  is  highly  oscillatory  and 
very  complex.  As  the  particle  grows  still  larger,  its  behavior  in  the 
radiation  field  approximates  that  of  a  large  solid  object  casting  a 
shadow  field  in  proportion  to  its  geometrical  cross-section.  The  Hie 
theory  includes  both  Rayleigh  scattering  for  snail  radii  and  the  **non- 
s elective"  scattering  in  the  limit  of  large  particle  radii.  The  theory 
is  covered  in  some  detail  in  the  next  section  cf  this  paper. 

Once  one  has  both  the  atmospheric  density  and  data  on  the  aerosol 
composition  and  size  distribution  as  functions  of  altitude,  location, 
and  time,  the  transmittance  can  be  computed  on  a  real-time  basis  by  some 
rather  sophisticated  computer  nodels. 

Aerosols  are  colloidal  particles  dispersed  and  suspended  or  falling 
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slowly  in  the  gaseous  mixture  cf  the  atmosphere*  Sone  examples  are 
smokes,  haze,  clouds,  small  droplet  (less  than  1y*n)  fogs,  and  finely 
divided  soils.  These  aerosols  have  as  many  sources  as  they  have  diff¬ 
erent  forms  and  compositions.  The  chemical  reaction  of  various  atmo¬ 
spheric  gasses  whose  minor  components  are  nitrous  oxides,  terpenes,  and 
hydrocarbons  can  produce  solid  particles.  The  presence  of  combustion 
products  from  inefficient  burning,  which  cay  be  locally  critical,  and, 
on  a  much  larger  scale,  from  forest  fires  and  volcanic  action  serves 
only  to  increase  the  atmospheric  load  of  particles  that  may  take  years 
to  settle.  Over  large  bodies  cf  water,  the  dispersion  of  solutions  into 
the  near  atmosphere  region  brings  about  the  formation  of  fogs  and  rain 
as  well  as  clouds.  As  a  class,  aerosols  may  either  help  or  hinder  the 
researcher.  These  particles  may  be  used  as  tracers  enhancing  a  received 
signal  to  detect  atmospheric  motio:  a,  distributions,  and  chemical  com¬ 
position.  However,  they  primarily  provide  a  large,  random  source  of 
environmental  noise.  In  this  latter  mode,  the  particle’s  presence  tends 
to  confuse  the  received  signal  at  the  sensor  and  provide  noise  that 
originates  outside  the  system  electronics,  making  it  difficult,  if  not 
impossible,  to  supress.  [6}  Historically,  aerosols  have  been  divided 
into  classes  to  *  facilitate  their  theoretical  modeling.  These  models 
involve  several  free  parameters  that  are  designed  to  “best  fit”  the 
observed  data.  The  classes  are: 

1.  Maritime  -  This  group  is  characterized  primarily  of  sea  salt 
with  some  water  of  hydration  and  a  chlorine  concentration  slightly  higher 
then  sea  water.  Salt  particle  concentration  decreases  rapidly  above  a 
height  of  about  500  m. 


UpWlfPPttpllMWHVWM  ■  •rnmnmmmtnmm . . . 


2.  Rural  -  The  largest  number  of  components  are  silicon,  iron 
and  some  sulfates  with  approximately  105  to  15*  organics.  Of  this  total, 
about  1/3  of  the  aerosol  i  ygroscopic,  whereby  air  moisture  is  readily 
absorbed.  This  same  class  is  also  found  at  sea  at  higher  altitudes  with 
a  large  particle  density  falling  off  rapidly  as  distance  at  sea 
increases. 

3.  Urban  -  Combustion  products  and  the  issue  from  industrial 
processes  make  up  the  bulk  of  this  class. 

4.  Troposphere  -  This  is  an  atmosphere  division  lying  above  the 
surface  boundary  layer  and  consists  of  the  rural/maritime  classes  with¬ 
out  large  particles. 

5.  Stratosphere  -  Primarily  {10%,  to  9C5)  consisting  of  the  large 
sulfate  ions  and  particulates,  this  group  also  shows  a  periodic  large 
influx  during  volcanic  activity. 

The  typical  atmospheric  analysis  technique  is  to  use  tvc-ended  sys¬ 
tems  requiring  a  source  of  radiation  of  known  properties,  a  propagating 
path  through  the  atmosphere  of  specific  length  and  reasonably  known 
conditions,  and  a  stable  receiver  at  the  path’s  end.  Experiments  such 
as  these  are  able  to  measure  aerosol  size  and  composition  distributions 
to  check  out  the  various  composition  and  scattering  models.  They  also 
predict  the  performance  characteristics  of  a  partial1  ar  optical  device. 
However,  two  well- separated  platforms  are  required  with  the  attending 
problems  of  phasing  and  communications.  Also,  only  low  altitude,  hori¬ 
zontal  work  can  be  done.  A  one-ended  system  would  eliminate  the  need 
for  one  of  th9  platforms  and  enable  a  full  altitude/azimuthal  capability 
in  a  real-time  environment.  There  is  a  system  that  can  easily  be 
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adapted  to  the  one-ended  process  and  is  available  today.  The  laser 
radar,  JT  ,  system  is  similar  in  operation  to  Radar  with  optical 
wavelenc  .3  instead  of  radio  waves.  In  fact,  the  acronym  is  somewhat 
the  same,  meaning  "Light  Detection  and  Ranging".  Conceptually,  a  laser 
pulse  is  emitted  into  the  medium  to  be  analyzed.  This  pulse  interacts 
with  the  components  of  the  medium  and  is  both  scattered  and  absorbed  in 
distinct  ways  depending  on  the  material  type  and  density.  Some  portion 
of  this  pulse  is  returned  as  an  echo  to  the  point  of  origin  where  a 
telescope-like  receiver  is  coupled  to  a  photodetector  that  converts  the 
received  optical  radiation  into  a  proportional  electrical  signal.  The 
pulse  return  timing  is  gated  to  the  range  of  the  atmospheric  cell  under¬ 
going  analysis.  The  intensity  of  the  echo,  the  echo’s  polarization 
change,  and  any  measured  frequency  shifts  all  undergo  involved  processing. 
The  outputs  can  consist  of  composition,  type,  and  size  distributions. 

Some  state  of  the  art  applications  of  the  technique  include:  raetero- 
logical  LIDAR  investigations;  middle  atmospheric  LIDAR  studies;  tropo¬ 
spheric  chemistry  and  diffusion  research;  studies  of  atmospheric 
propagation  and  radiative  transfer;  absorption,  Raman  and  fluorescent 
spectroscopic  applications;  and  pressure-temperature  measurements  of 
atmospheric  layers.  C23  A  recent  theoretical  paper  indicated  that  a 
series  of  measurements  involving  optical  scattering  and  absorption  can 
yield  about  7  or  8  independent  pieces  of  information  concerning  the 
aerosol  size  distribution  within  the  analyzed  cell,  covering  a  range  of 
particles  with  radii  from  O.l^m  to  5/<m.  £2q| 

Consider  such  a  single-ended  LIDAR  analysis  system  on  board  the  type 
of  mobile  platform  in  use  today  by  both  scientific  researchers  and  the 
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II.  GENERAL  SCATTERING  THEORY 


Scattering  is  the  process  whereby  an  incident  electromagnetic  wave 
loses  energy  and  some  portion  of  this  energy  undergoes  re-emission.  In 
general,  the  scattered  wavefront  has  a  different  propagation  direction 
than  the  incident  wave.  See  figure  2. 

For  particles  with  a  low  manner  of  electrons,  the  scattering  can  be 
computed  using  quantum  mechanics.  In  this  sense,  the  incident  wave  can 
only  interact  with  the  stationary  states  oi  the  scattering  particle's 
electronic  system  within  a  small  bandwidth  of  frequencies  about  tho 
allowed  transition  levels.  One  can  discuss  Rayleigh  scattering  as  the 
absorption  of  the  correct  energy  photon  and  the  subsequent  re-emission 
of  a  photon  of  essentially  the  sane  energy.  An  energy  level  diagram 
uould  look  like  figure  3. 

Likewise,  Raman  scattering  involves  a  definite  frequency  shift  be¬ 
tween  the  incident  and  scattered  radiation.  If  the  scattered  frequency 
is  less  than  that  of  the  incident,  the  process  produces  the  Stokes 
lines  spectroscopically.  And  when  the  reverse  is  true,  the  anti-Stokes 
condition  exists.  DO  See  figure  4. 

However,  for  molecular  and  larger  particle  scattering,  the  collection 
of  electronic  states  becomes  great  enough  to  let  us  deal  with  the  prob¬ 
lem  classically.  The  particle  is  considered  to  be  a  solid  dielectric 
with  a  complex  index  of  refraction. 

One  of  the  most  widely  studied  problems  in  diffraction  theory  is 
electromagnetic  and  acoustic  wave  scattering  by  a  sphere.  The  problem 
of  diffraction  of  a  plane-monochromatic  wave  by  a  homogeneous  sphere 
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characterized,  by  its  radius  and  complex  index  of  refraction  was ’ rigor¬ 
ously  solved  by  G.  Mie  in  1908.  His  analysis  involved  the  solving  of 
Maxwell’s  equations  for  the  electromagnetic  field  when  an  incident  plane 
wave  interacts  with  a  material  interface.  When  the  discrete  boundary 
separates  media  of  different  optical  characteristics,  a  scattered  wave 
is  generated. 

Today,  an  extensive  bibliography  of  approaches  to  this  problem  exists 
The  author  will  follow  the  analysis  presented  by  Kerker  (1969)  03 
throughout  this  section  in  order  to  present  a  consistent  derivation. 
First,  the  salient  features  of  the  Maxwell  theory  wall  be  listed,  lead¬ 
ing  to  the  wave  equation.  The  solution  satisfying  the  appropriate 
boundary  conditions  will  be  obtained;  the  resulting  inf  ini-'  e  series 
solution  not  only  completely  defines  the  scattered  wave,  but  gives  the 
electromagnetic  states  of  the  particules’  interior  as  well.  The  single 
particle  scattering  results  will  be  compiled  and  extended  to  cover 
collective  scattering. 

Maxwell’s  equations,  in  rationalized  MKS  units:  D2L03 
Coulomb’ s  law  •  Q  s  <°  ( 1 ) 

Ampere’s  law  TxE  ^ 

Faraday’s  law  y  x  E  (3) 

absence  of  free  magnetic  poles  V  sq  (4) 

The  symbol  jl  denotes  a  vector  quantity.  Following  are  definitions 
of  vectors  used  in  derivation: 

Z  -  electric  field  intensity  (volts/meter) 

U  -  magnetic  field  intensity  (anps/meter) 

jD  -  electric  flux  density  (dielectric  displacement,  coulomb/meter 


& 


is 


1, -  magnetic  flux  density  (magnetic  induction,  weber/meter  ) 

2 

J  -  electric  current  density  (amps/ meter  ) 

'X 

/°  -  electric  charge  density  ( coulomb/ meterv } 
t  -  tine  (seconds) 

These  equations  './ill  define  the  electromagnetic  field  configuration 
for  all  points  in  space  and  uniquely  determine  the  field  vectors  given 
the  media  characteristics  below. 

E-'SE  (5) 

'~'*L  (6) 

Z  -  S'  &  (7) 


The  media  parameters  thus  defined  are: 

£.  -  permittivity  (electric  inductive  capacity,  farads/meter) 

-  permeability  (magnetic  inductive  capacity,  henrys/neter) 

6  -  conductivity  (specific  conductance,  l/ohn-meter) 

In  free  space  (vacuo): 

£»  3.8542  X  10~12  farad/m 

n 

/+*/*>*  10“  henry/n 

As  a  consequence  of  the  invariance  of  the  Maxwell  equations,  the  quantity 
C*er  a  has  the  units  of  velocity  and,  again,  in  free  space: 

C»  fatf*  »  i.W1/5  X  |0*  m/sec  (S) 

This  last  quantity  is  the  speed  of  light  in  vacuum.  ’.Jhen  the  quantities 
andefare  independent  of  direction,  the  region  in  space  they 
characterise  is  said  to  be  isotropic.  Two  more  equations  are  required 
to  complete  the  picture. 

Lorentz  Force  £  r  +  X  ®)  (9) 

lievton’ s  second  la u  of  notion  £  -  ma  (10) 


'ihere:  ?  -  force  on  charged  particle  (newtons) 
ju  -  velocity  of  charged  particle  (a/s) 
jn -  mass  of  charged  particle  (kg) 

&  -  acceleration  of  charged  particle  (a/s  ) 

All  macroscopic,  linear,  electromagnetic  phenomenon  are  described  by 
equations  (l)  through  (10)  with  the  associated  quantity  definitions. 

The  classical  dynamics  of  interacting  charged  particles  and  electro¬ 
magnetic  fields  are  thus  specified. 

A  certain  degree  of  convenience  is  obtained  using  the  definitions  of 
the  two  following  auxiliary  vectors: 

P  =  2-f.I  (11! 

02) 

2 

Where:  ?  -  electric  polarisation  (c/m  ) 

M,  -  magnetic  polarisation  (A/r.) 

One  can  see  that,  in  free  space,  these  quantities  vanish  identically  and 
they  may  then  be  considered  as  a  measure  of  the  effect  of  matter  on  the 
local  fields.  Some  dimensionless  quantities  may  also  be  defined: 

Kc  =  §/e0  03) 

Km  =  ^Ao  04) 

Where:  f(C  -  specific  6  (relative  dielectric  constant) 

-  specific^ (relative  magnetic  permeability) 

The  media  may  further  be  described  v/ith: 

£  -  os) 

m  =  x„h,  (i« 

’■/here:  electric  susceptability 

magnetic  susceptability 

The  orevious  set  of  eauations  can  be  combined  to  obtain: 
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Xe=  Ke  -1 
X m~  K*m  "i 


(17) 

m  (1S) 

Maxwell’s  equations  apply  to  regions  where  the  media  may  be  repre¬ 
sented  by  the  continuous  functions  £,/*,  and  O'.  1/hen  a  medium  discon¬ 
tinuity  is  encountered,  such  as  the  conditions  across  a  regional 
interface,  the  fields  must  satisfy  a  specific  set  of  relations  at  the 
boundary.  See  figure  5.  Let: 

£5  -  surface  current  density  on  interface  S  (A/a^) 

2 

/£  -  surface  charge  density  on  interface  S  (c/n  ) 

The  Maxwell  equations  (l)  -  (4)  are  cast  in  point  form  for  each 
differential  point  in  space.  One  can  apply  various  integral  theorems 
to  change  the  same  equations  to  integral  form  and  use  them  to  define  the 
conditions  that  the  fields  must  satisfy  by  crossing  the  interface,  jj7] 
The  boundary  conditions  obtained  are  inherent  in  the  Maxwell  equations. 

Boundary  conditions: 

'  the  tangential  component  of  is  continuous 

fe  -  £|)  x  2  =  O  (19) 

(b)  The  tangential  component  of  H  undergoes  a  discontinuity  equal  t_ 


0U-«i)*  o=  Is 


(20) 


(c)  The  normal  component  of  Jj,  undergoes  a  discontinuity  equal  to  /& 

'2,) 

(d)  the  normal  component  of  is  continuous 

CS*-S.)tt  =  o  <22> 

For  all  cases  of  dielectric  scattering,  it  is  assumed  that  there  is  no 
free  surface  charge  or  current  density.  This  assumption  begins  to  break 
down  for  large  particles  in  the  presence  of  storm  activity.  However,  we 
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will,  for  simplicity,  assume  that  =  0,  and/°s  s  0.  Therefore,  the 
tangential  components  of  and  H  and  the  normal  components  of  D  and  3^ 
are  continuous  across  the  boundary. 

Charge  conservation  is  included  in  the  Maxwell  theory  and  is  written 
as: 


continuity  equation 


V  '  +  If  /?  s  o 


(23) 


O 

Where:  -  volume  current  density  (A/nr ) 

-  volume  charge  density  (c/m  ) 

Throughout  the  following  derivation,  a  few  vector  identities  will  be 
required. 

Note:  A  -  continuously  differentiable  vector  field 
0-  continuously  differentiable  scalar  field 

r*(vxA)=  o 

V  =  O 
=  a  .  (V0)  +  p'(v-d) 

VxOar>l)  ^(vx4)  +  (^)xA 

V  X  (rx4)  =  V(V-£)  -  V2A, 

If  the  fields  described  by  (l)  -  (4)  are  invariant  under  a  gauge  trans¬ 
formation,  then  they  are  well  defined  by  the  preceding  formalism.  pa 
We  define: 

A  -  vector  notential  field 

/v* 

0  -  scalar  potential  field 
in  the  following  relations. 

B  =  VXA,  (25) 

Compare  this  to  equations  (4)  end  (24).  Using  equation  (25)  in  (3),  we 
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get  the  following: 

V'iZ  =  “Vx(|fA) 

=>*  X  (e  +jt  =  0 

How,  refering  to  (24), 

=>  £  +  ir  A  » 


nz 

£* 


,thus 


(26) 


A  and  0  are  not  unique.  Let  g  be  any  scalar  function,  then  the  new 
potential  can  be  defined  as: 

=  ^  -VQ  (27) 

(23) 

Now  the  fields  2  and  3,  are  reconouted. 

A»  mJ  * 

31  =  F  X  (A-F9)  =  £  -TX  (F9^  *  g 
=>  O'  =  8 


e  =  £ 


The  physical  observables  are  the  fields  2  and  3_ and  it  can  be  seen 
that  the  gauge  transformation  defined  by  (27)  and  (23)  leaves  these  fields 
unchanged.  Thus,  the  Maxwell  equations  are  gauge  invariant. 

We  nov  turn  fron  formalism  to  a  more  practical  aspect  of  electro- 
nagnetic  theory  and  introduce  the  Wave  Equation: 

rxd):  vx(rx£)-  ^  (y  x  s’) 

Note  that  the  tine  and  space  operations  of  differentiation  are  connutative 
for  these  continuously  differentiable  functions. 

OO, (2*0:  v  (y-E)~ 

The  assunption  is  always  Bade  that  all  space  and  tine  derivations  of  the 
paraneters  £,/*,  and da re  either  slowly  varying  with  respect  to  the  rate 


24 


of  change  of  the  fields  themselves  cr  identically  equal  to  zero. 

=>  v(v£)-'V1£*  ~/*4t  x  ^ 

(0,  Cz),  (5),  (7)  ;  r  •  (&  £  )  =  t  v  •  E  =  P 

=>  v  •  £  =  /% 

£je£  + 

->  v  x  £  =  o'  £  +  ^It  & 

->  vC/Vfi)-  V2g=  v^lh  €jt 

The  sane  argunent  also  applies  to  the  quantity  fij £  concerning  its  rate 
of  change  with  the  rate  of  change  of  the  fields. 


=>  V  E. 


Therefore:  V 


-  f  E 

dt  **  3t*- 

*£  _  r/<iE  -  t/.ai  E  =.  O 

^  3t  A/ 


— —  ,  —  ^  ~  '  3d1"  ^ 

This  equation,  (29),  is  the  well-known  nonhoaogeneous  equation  for  danced 

wave  notion  -  the  lave  Equation.  Similarly, V  X  (2): 

vx  (vx  H)  =tt  Ox &•)+-(.  TXJ) 

(5),  (i),( 2*1):  v(v-jjVv1rii£-fe-Oxl)<.^TX£-) 

(H)  — >  V  •  (>**■  H  1  V-  ti  n  0=r>  V«|4  =  0 
(3):  -VZH  = 

Cfc):  = 

— >  V  hi  -  1  =  O  (30) 

The  relations  for  non-dissipative  media  are  obtained  by  setting  cf  -  0. 
3oth  equations  (29)  and  (30)  actually  represent  3  equations  each.  Each 
component  of  the  vector  field  must  satisfy  these  equations  for  the  vector 
field  itself  to  do  so.  Let  the  symbol  d  correspond  to  any  of  the  six 
components  of  the  two  fields  j3  and  ij.  The  scalar  wave  equation  can  be 
expressed  by: 

VJU-cr/uia._eA4^  =0  (3,) 
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The  fields  described  by  (31 )  represent  a  transfer  of  energy  in  both  space  .  ; 

and  tine.  The  electromagnetic  wave  has  a  tine  periodicity  that  sug-  ; 

gests  a  forn  for  the  energy  flux  crossing  a  unit  area  per  unit  time: 

S  *  £  *  H  (32) 

Equation  (32)  is  the  Poynting  vector  equation,  whose  tine  averaged 
magnitude  over  the  period  of  the  wave  yields  the  wave  intensity, 

1=  '/rjjildt  os) 

where f is  the  period  of  the  wave.  The  field  symbols  E  and  H  in  (32)  are  * 

the  real  parts  of  the  complex  fields  whose  components  satisfy  (31). 

The  form  of  the  scalar  wave  equation  suggests  a  harmonic  time  de- 

I 

pendence  of  the  form:  * 

U  (c,t)  =  -PCc)elU,t  (34)  i 

I 

where  co  is  the  angular  frequency,  i  =VM),  and  r  is  the  usual  three- 

} 

dimensional  position  vector.  Hote  that:  t 

2^"  =.  lcoU.  and  '%'^i  -  U 

then,  (31):  r*U-«^c»U  =  0 

I 

=>  VZa  +  Z  U)z  -  i/s  cr  Co)  u  =  O  i 

Define  the  propagation  constant  1c  by:  \ 

K2  (35) 


then:  V  2U  ¥  RZ  U  -  0 

(36) 

Equation  (36)  is  the  Helmholtz  Wave  Equation. 

Another  form  of  (35)  is 

K  a  —  l  43, 

(37) 

where:  ,  =  to  f  ^[(l  H}  ^ 

(38) 

*•  -- 

(39) 

It  is  the/3,  term  that  determines  the  damping  undergone  by  the  wave  per 
unit  path  length;  the  energy  is  dissipated  as  Joule  heat  in  the  medium 
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•.men  C  >  0 


Recalling  that  the  quantity 

v  *  (/*£ )/z 


(40) 


is  a  velocity  which,  in  free  space  is  c,  the  relation  for  an  electro¬ 
magnetic  wave  in  a  medium  described  by  light  speed  V  is 

*  V  (41) 

where:  X  -  wavelength 

h)  -  frequency 

V  -  velocity  of  the  wave 

Also, . 2Tl'v  :  u)  '  ...  (42) 

and  this  is  the  angular  frequency,  u> ,  mentioned  earlier.  In  free  space, 

(43) 

The  index  of  refraction,  n,  is  defined  as 

n  =  c/y  (44) 

X  x  / 

U1),  (43),  (44):  7A 

=>  r\  A  =  A0  (45) 

It  is  clear  that  A 0  is  the  wavelength  when  the  wave  is  prop -gated  in 
free  space.  In  the  same  medium,  O'  -  0  implies  that 

~  0  =>  °f  I  S.  UJ  (/*fl  £0)  »  o)/c. 


Ae^c 


Koiu)/c=  E1Ti)/a0V 
=>  Ko=  2ir/^o 


(46) 


This  kQ  is  the  free  space  propagation  constant. 


The  complex  refractive  index,  m,  can  be  defined  as 

m  =  k/k0  =  n  0  -  i  K)  (47) 

’..'here:  X  ~  index  of  absorption  (extinction) 

As  an  example,  suppose  that  there  is  a  plane  wave  propagating  along  the 


z  -  axis  whose  x  -  coraoonent  is: 

£,~a. 


(47 

=>£„./!.  “*'**•%-"«.**■ 

■=*>Ey  = 


-  e 


^  -inK0  (/-  L  X)  z 


-tlk-)C2 


The  attenuation  of  the  wave  is  seen  explicitly  as  e  “  ®  .  An  expansion 

of  this  example  shows  that  the  resultant  plane  wave  solution  leads 
directly  to  the  Stokes  parameters.  Assume  the  propagation  direction  is 
zf  as  stated  earlier,  and  that  the  monochromatic  wave  is  traveling 
unbound  in  a  source  free  isotropic  medium.  It  is  assumed  that  the  wave 
can  by  completely  specified  in  the  xy-plane  by  the  z, t  coordinates  of 
the  wave.  Thus: 

r2u  =  ily.  (*.*)  =>  lu  Jut 

b  X  “  5  V  -  ° 

(34),  (36):  aiut  =Q 

d-  V  .in 


l  + 


k*T  * 


oz 

Substituting  a  trial  solution,  f  =  ae"  , 

=>  ab*eb%  Kltt*b*-s,o  =?  t"-=  -k%  b*  1  ik 

=.>  f  (z)  =  ae  i  iKz 

From  the  definition  of  k,  -iks  is  selected  as  the  value. 

.  (v  xh  -  iCoit-Kz) 

U  (Z,W  -  CLZ 

This  equation  describe  the  form  for  each  of  the  six  vector  components 
of  S  and  H.  Recalling  (l),  under  the  present  constraints: 

(,)  f  MZ  =  o  mi  U= 

^t,  la.  du  ^  c.  -  o 

'0l 

=>e,=  4«‘^  5  •**  £h-  4 e 


(3),  (6) 


vxfs  --£iit  ->  -sj*Mx  -  Z^-Y  i 
at  x  ”  3z  y  ’ 


-<£&£»  Alt  ~d±M  z  -n  __s  (4  r  n 

=>  df  y'<5Z*  >  c)t  u  _>  Hz  0 

It  is  quite  sufficient  to  describe  the  behavior  of  J2  as  the  value  of  I{ 

can  be  readily  obtained  using  equations  (l)  through  (4).  Given  the  fora 
of  u(s,t) : 


3 

SH 

3i 

Y 

-  i  ui  Hy 

11 

1 

r*» . 

X: 

) 

<5  £ 

a  z 

■y 

—  - 1 

n£i 

+  i  k.  £y 

=  ? 

— 1C 

£a 

"  l/c  £  % 

=3> 

Hi 

: 

ft 

yCS  Ui 

rvAx  =/"wA 

reduced 

using  (35) 

to 

the  following 

=>  ^  (f^r) * 

Defining  the  intrinsic  impedance  of  the  medium  for  plane  waves  as 

z°  =  (^%)  4  (43) 


rx  l£y 


H  y  ”  «x  -  Z° 

Finally,  ^  for  the  plane  wave  is: 

-*z  +S,) 


E*  =  Ae 


£w  =  r  (49) 

£*  =  0  J 

where  and  a^are  arbitrary  phase  angles  to  account  for  the  general  his¬ 
tory  of  the  wave.  The  H  wave  is  also  determined  by  eauation  (49): 

[  (5=) 


nz  =  O 

lie  see  that: 


E-H  =M*-EvHy  -  +■  ZipL  -o 

^  ^  O  ^  ft 


H 
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(51) 


=?>  E  •  H  :0 

Thus,  and  v  are  perpendicular  to  each  other  for  the  plane  wave 
solution. 

The  concept  of  polarised  light  can  now  be  introduced  with  the  help 
cf  equation  (49), X, ,  and^.  The  general  plane  wave  solution  describes 
an  elliptical  rotation  of  the  2  vector  in  the  olane  perpendicular  to  the 

M  * 

direction  of  propagation.  Defining  the  phase  difference  &sJ" : 

=  (52) 

The  representation  of  the  state  of  a  polarized  wave  can  then  be  uniquely 
determined  using  the  Stakes  parameters,  which  are: 

5o  =  A2+  B2  =  A,2  +•  0?  (53) 

S,s  A2-B2  =  S0  COS  2  V*  COS  2  x  (54) 

s2=  2 C05  /  =  30  S»M  2  V'  Cos  2  x  (55) 

s3=  2AS  SINcT  s  S0  SIN  2  X  (56) 

The  lengths  and  31  and  the  angles  ar.dX  are  illustrated  in  figure  6. 

The  quantity  (mt  -  kz)  can  be  eliminated  from  (49)  using  equation 

(52),  where  E  and  E  refer  to  the  real  marts  of  the  conolex  wave: 

f  [  5f)  "  2£*%  COS  5  =  SIN2  S  (57) 

VA/  \B/  As 

This  is  the  equation  of  an  ellipse.  The  ellipse  in  figure  6  is  drawn 
"right-handed",  that  is,  the  head  of  the  electric  field  vector  follows 
the  arrows  when  viewing  the  wavefront  face  on.  In  order  to  exhibit  right- 
handed  polarization,  sin$>0  and  0<X£ff/h  Linear  polarization  exists 
when: 

§  jaO.HtZ.iS,..*  (56) 

and  circular  polarization  corresponds  to 

§  =  j  =  £l,i3,is,t7J  ...  (59) 
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Pictorial  definition  of  the 
quantities  ,  3. ,  'f  ,  X 


2A 

(0  6  *<ir) 


Figure  6 


The  Stokes  parameters  satxsfy: 

S0  x  5,  +  S2  +  53  (60) 

and  a  complete  state  of  polarization  is  obtained  from  only  three  in¬ 
dependent  parameters. 

Natural  light  is  defined  as  the  state  where  the  intensity  of  light 
is  the  sane  in  any  direction  perpendicular  to  the  propagation  vector,  k. 
The  phases  of  the  various  components  in  natural  light  are  randomly 
varying  so  that  no  correlation  exists  between  them. 

$?  *  =  S3  =  0  (61) 

The  combination  of  elliptically  polarized  beams  can  be  accomplished  using: 

S0  =  .<f  Vc  J  S,»  i  s,l  ,  S2i  £  5 21  }  $3  =  £  s3l  (62) 

Then  for  non- identical  polarization  states  for  all  beams: 


S*  >  s?  +  si  *  si 


(63) 


Partially  polarized  light  can  be  described  as  a  superposition  of  natural 
light  and  a  set  of  incoherent  beams: 


Z  >  S3  *  S3  3 


and 


(64) 


S0=  Sft  +  So  ,  S,  S 

S t  ■-  (sl  +  s^sj)* 

Finally,  the  fraction  of  polarization,  F,  can  be  defined  as  the  intensity 
ratio  of  the  polarized  part  of  the  wave  to  the  total  wave  superposition: 

F  zZ^t  s  So/CsJ  +  S t)  Co^F^O  (65) 

It  can  be  shown  that  the  parameters  of  a  scattered  wave  are  linearly 
dependent  on  the  parameters  of  the  incident  wave.  03  This  relation 


is  written  as  follows: 

- 


is  =  o/RMfii 


(66) 


where:  R  -  observation  distance  to  scatterer 
i  -  incident  wave 


s  -  scattered  wave 


The  modified  Stokes  parameters  are  defined  first  by  writing  equation 


(49)  j 

as: 

iuit 

LUit 

’  £v 

(67) 

then, 

I.  -  K 

I2 

(63) 

X2  -  |  £j 

J2 

(69) 

U=  2  R 

(7C) 

V  -  2  I 

(71) 

The  444  Stokes  matrix  ej,  looks  like 

this: 

if„r 

If,,!2 

G  - 

iy2 

ifj2 

(P2|1jD 

(72) 

*  -WfAS-p.jyO 

2  Xm('Pitfji) 

S- 

2  ImfriifaJ+Wz 

i)  Re 

where  the  f. .  are  defined  in  terms  ox'  the  scattering  functions  given  in 
m-  For  spherical  particles,  the  quantities  =*  f21  ~  0  50(1  ® 
contains  four  constants. 

Having  laid  the  foundation  for  the  Hie  treatment  of  scattering 
phenomenon,  let  us  now  consider  a  plane  electromagnetic  wave  incident  on 
an  isotropic,  homogeneous  sphere  of  arbitrary  size.  The  object  interferes 
with  the  free  propagation  of  the  wave  and  secondary  waves  are  generated 
that  are  predictable  by  Maxwell* s  equations  and  the  boundary  conditions; 
and  these  secondary  waves  form  a  scattered  radiation  field.  The  incident 
wave  should  be  linearly  polarized  for  two  reasons.  First,  if  the  light 
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is  unpolarized,  then  it  can  be  resolved  into  two  linearly  polarised 
components  each  of  which  act  independently  of  the  other.  One  only  needs 
to  consider  the  interaction  of  the  sphere  with  a  plane  (linear)  polarized 
wave  and  then  apply  the  principle  of  superposition  to  obtain  the  total 
effect.  The  LIDAR  laser  pulse  is  plane  polarized,  by  definition,  and 
this  is  the  wave-type  that  is  considered  for  tropospheric  probing.  A 
geometric  diagram  of  the  problem  is  shown  in  figure  7,  suggesting  a 
spherical  coordinate  system  as  the  most  convenient  for  analysis. 

The  fields  of  electromagnetic  energy  in  the  region  of  space  filled 
by  the  particle  and  the  incident  wave  are  resolved  into  three  components: 

a.  the  incident  wavs  -  S.,  H. 

b.  the  Darticle  interior  wave  -  E  ,  H 

*  -r  ~r 

c.  the  scattered  wave  -  E  ,  H 

**  S 

For  these  quantities  to  properly  define  the  fields,  they  must  each 
satisfy  Maxwell’s  equations  (1)  -  (4),  the  Have  equations  (29),  (30), 
and  the  boundary  conditions  (19)  -  (22).  The  internal  field  S  .  H 
must  match  the  external  field  E.  +  E  ,  H.  +  H  in  accordance  with 
these  relations. 

Kerker{]2)et  al.,  suggests  that  the  derivation  deal  with  the  scalar 
wave  equation,  (31 ).  Two  new  functions  may  be  introduced,  the  electric 
Hertz  vector, IT] ,  and  the  magnetic  Hertz  vector, |£,  and  are  defined  as 


Bi  =  /£  v  U» 

(73) 

E,=  r(v ■%)-/<£&  u, 

(74) 

D2=-^£V  x^-TT,  ^ 

(75) 

Hj=  V(v-Tri)-S*z£p-I]'z 

(76) 

follows: 


Suitable  transformations  may  be  applied  tc  define  other  Hertz  vectors 
as  long  as  the  fields  are  left  invariant.  The  vectors  thus  defined  must 
satisfy  the  vector  wave  equation: 


where  the£  andji  result  from  distributions  of  electric  and  magnetic 
dipoles.  Therefore,  two  vectors  can  be  used  to  describe  all  the  fields 
present. 

The  transverse  magnetic  wave  (Til),  also  called  the  electric  wave, 
has  a  zero  value  for  the  radial  component  of  its  magnetic  field  intensity, 

-  0.  Likewise,  the  transverse  electric  wave  (TE),  the  magnetic  wave, 
has  a  zero  radial  component  of  electric  field  intensity,  E^  -  0.  The 
TM  ’wave  is  conceived  as  being  the  result  of  oscillating  electric  di¬ 
poles  in  a  spatial  region,  and  the  TE  wave  as  the  result  of  oscillating 
magnetic  dipoles.  The  total  solution  is  a  construct  of  these  two 
wavefronts. 

The  Hertz-Debye  potentials  are: 

“  V  *  1?  ~  %  (79) 

-  v  •  Trz  =  ir2  (so) 

These  potentials  are  the  solutions  of  the  scalar  wave  equation.  Solving 
(79)  and  (SO)  for  the  potentials  and  applying  the  appropriate  boundary 
conditions,  one  can  then  find  the  field  vectors  describing  the  TO  and 
TE  waves.  The  total  field  may  be  derived  by  the  direct  addition  of  the 
component  waves. 

The  following  equations  for  the  total  field  vec” .  components  can  be 
obtained  by  applying  these  concepts  to  the  above  relations.  In  spherical 
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coordinates,  we  have: 

Er>E,r*£2r  =  (rir,)+  KlCrl0*O  (S1) 

©  -  T  Cr,,i)  +  K2>  Sift©  Cr,0  (-82) 

+  ^-rsift©  (33) 

-  Hir  +  Hzr  —  O  '*'^1  (r  1U^+  K1  (34) 

H©s^i©  +  H2q~  ram©  ^)  + r  3r^©  ("r  '0  (35) 

+  K,  -p  jg(rii,V"r Jr 3^ C*  ’*' 2^  (36; 

with  the  propagation  constant  defined  in  the  following  fora  by  (35): 

Ic2  =  “K,  Kz  (87),  where  K  ,=  LuJt  +  o’  (?s),  \CZ  =  iu>  (35) 

^  has  been  dropped  as  all  redia  are  assured  to  be  non-nagnetic. 

Recalling  the  scalar  wave  equation  for  a  function  with  haracnic 
tine  dependence,  e~  °,  we  have: 

(36)  r2u  +  K2U  =  o 

(34)  U=f(r)e"‘ 

The  spatial  part  f  then  satisfies: 

t2f(r)+K^(d:0  (90) 

The  llertz-lebye  potentials  are  constructed  such  that  they  are  solutions 
to  this  wave  equation.  LetTT'(r,  0,  <f>  )  -  f(r)  represent  cither 
potential,  7T,  or 7£z 

v2ir  +  k’-it  =  o  (9i) 

This  equation  is  the  one  that  nust  be  solved;  in  spherical  coordinates: 

r  (“tr*  rir)  +  r*slne  ii  (SLneje)+  r\s.n\*  K  V  s  O  (92) 

It  can  be  solved  using  the  technique  of  separation  of  variables,  where 
ITis  assured  to  be  a  product  function: 

1TCO  =  RCr)®  (©)§(«*)  (93) 

s>  T7*  -  ®  5  ^  r  ^  (94) 


I 

ii 


i  -I 
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W: 

J"  i 


/ri-\ 

\">  f 


(96) 


=>  &(smeI)  =  R2f.(s,"e!§) 

=>^rt=  RO$c 

(92)  be=o=eS:®J  £0  +7fi (sm  a  =  -K1R®1 

3ir.cs  IT  ^  3,  then  *:c  multiply  by  r2/ir: 

x  ' —  i./ siM©  a®)  J— -  ili  _  _k*  2 

R  3 rl  ©sins  3©y+$siM*o  3  jsf  4  -  *•  r 

v  r  i!  r  FU  Klr 1  Q - X.  /sin  e  365)  J _ (=■?) 

R  3rx  -  "  {©swe  a©\  ^tIsim10 

This  equation  is  equal  to  a  constant  which  ray  be  checked  by  talcing  the 
derivatives  of  both  sides  with  respect  to  r,  ©  ,  <P  .  This  constant  has 
the  forr  of  r.(n  +•  1): 

^.^rRVKV2  =  n(n*0 

Thus  the  radial  ecu at ion  is: 


hi  rRCr-)  SkfPJ  rR(ri  =  O 


(93) 


"sing  (93),  ecuaticn  (97)  becomes: 

_J _ Jl/si^6  3®n  j _  ,3*5  ,  v 

©SIN©  38  (  "5©/  + $51**0  XSF*  s  “n^n1'0 

pSI*©  _d/$IN©  3©\  /  _JL-d*£ 

LST^  d©(  •5§)  +  n(n+0s.fl  ©J=  ^ 

By  the  sane  differentiation  argument,  each  side  of  this  equation  is 

2 

ecual  to  a  constant,  m  : 

‘  =>  *  «* 

Thus,  the  0 -angle  enuation  is:  __ 

fc*  fc(*eif£9+i"<"*0-  0 

Using  (98)  and  (99),  (97)  takes  the  form  of: 

»*•£  _  2 
2  T**  -  m 

Thus,  the  <t> -angular  ecuation  is: 

In  the  above  equations,  n  nay  take  on  any  integer  value  and  n  nay  be 
any  one  of  the  following  set:  -r.,  ...,  0,  n.  Zcuatior.s  (99)  and 


(99) 


(  1 0TI  ^ 
v « uo ; 


il 


! 'tililtul'i'llllMill lill'll !<< ij* ii<tii<;Hji'l) h'llill^illill]1  niliiljiinliiil’i 'p<<il ' i  111  l<i  iP H]n |i)Ji ,  ,]]ji  ilM'll'W  . . . . & . a . it . a . . . . & m . . to a . tarn . a . am . . . . . miMM . ti^ 


( 1  CO  ^  con  bine  to  form  the  Snierical  '’monies  which. are  defined : 

C  (0=C®>)fmW  00,) 

’..’here  A  is  a  normalisation  constant  to  ensure  orthonormality  of  the 
n 


The  solution  to  (700)  is  a  straightforward  combination  of  sines  arc 


*  2  . 

cosines  as  n  is  always  positive: 


J  (0)  -  dm  C os  O$0  +  jom  5in 


(102) 


legation  (99/  is  the  familiar  Associative 


one  rsuuiar  soiusicn: 


®Ce)  =  P"  ;  (cos  ©) 


{ I  p. « '* 

V  * 


..  /  /-s.  \  •-  ....  ,  ..  .  . 

•Mere  zze  '  (cos  &  }  are  wie  associative  legeriore  ceiyr.ctnaiS. 

Using  Rodrigues’  formula  (l]  : 

z'  n  iV1/,  ,^**1  Jrntn 

rn  (COS©)  —  ( 2  0  i  /  (l-X  )  -r  TST%  (104; 

where  :t  =  cos  ©and  the  negative  values  of  n  are  included,  -n  £  m  4  n. 
Uher.  a  -  C,  then  (9?)  is  the  leger.cre  equation  whose  editions  are  the 


legendre  polynomials: 


>„&>'.  (tV)-’  (£)  (x‘-')n 

Rcruation  (104)  may  be  rewritten  as: 

(I 


005) 


(106) 


file  complete  spherical  harmonic  solution  includes  the  normalisation 
constant  and  rewrites  the  azimuthal  solution  as  a  coupler  exponential.. 
These  functions  are  completely  orthonomal  as  well.  Thus: 

Y>,  <f)  £-  0"7w  |E§<  Pn(m\c«  «) « iw!r  (,c?) 

The  (-1 )n  term  is  a  phase  factor  defined  as  the  Condon-3 hortly  phase,  and 
nay,  in  some  texts,  be  dropped. 

(Rote:  in  the  remainder  of  this  paper,  X  will  use  j  instead  of  seta. 


mr m 


which  is  used  in  the  literature.) 


The  Radial  equation,  (93),  has  solutions  defined  in  terns  of  the 
Bessel  functions: 


■Xn(Kr)  =  -'/f  Kr  A/,«  Vj.  0<  rO 
=>  rRn(r)  *Cnl^(l‘r)<-dnK1.(Kr) 


(108) 

(109) 

(no) 


(111) 

(112) 

(113) 


These  are  the  Ricatti-3essel  functions,  where  J  , «  and  II  ,  i  are  the  half 

7  n+2  n+o 

integral  order  3essel  and  Neumann  functions.  By  using  the  infinite  series 
definitions  of  the  Bessel  functions,  one  can  '.cite  the  Ricatti- Besses 
function  in  the  following  derivative  fora: 

-(-Dl (ikrfzf*) 

where  z  -  x  +  iy 

The  function: 

=  f„  (Kr)  .  i  .■=  JTZ  (*,)  OU) 

vanishes  in  the  limit  as  Icr  approaches  infinity,  and  will  be  useful  in 
later  work.  The  function  H^N(kr)  is  the  half  integral  order  Hank  el 
function  of  the  second  kind. 

Now  it  can  be  seen  that  the  solution  to  the  scalar  wave  equation  is 
characterized  by  the  numbers  n, ‘m: 

rC°=  rRnM®<"%) 

and  (102),  (103),  and  (1 70) : 

=>  fC„WteUnCK<i>H?^  („5) 

The  general  solution  may  now  be  obtained  as  a  linear  superposition 
of  these  particular  solutions. 

rirco  =  U  rir£"V.e,<0  (n6) 


me  m-i-n 


Equation  (116)  is  the  general  solution  of  the  scalar  wave  equation  in 
spherical  coordinates. 


The  fields  in  and  about  the  particle  can  be  described  in  terns  of 
their  individual  pairs  of  potential  functions; 

a.  the  incident  wave  -  17/  3  TT/ 

b.  the  interior  wave  - 

s  $ 

c.  the  scattered  wave  -  irr>TP* 

The  sphere,  whose  origin  coincides  with  that  of  the  spherical  coordinate 
system  and  is  of  radius  a,  is  isotropic  and  homogeneous,  optically 
described  by  a  complex  propagation  constant,  : 

(47)  ~  wi  ^ 0  (117) 

The  isotropic,  homogeneous  medium  is  a  dielectric  described  by  a  real 
constant,  k^. 

Kz  -  (113) 

The  relative  refractive  index  is  defined: 

m=  K,/kz  *  rv\Jmt  Oi9) 

The  incident  plane  wave  propagates  along  the  positive  z-axis  and  is  con¬ 
sidered  to  be  of  unit  amplitude  with  polarization  axis  parallel  to  the 
x-axis. 

\h\'  |e'tk*Z|  020) 

In  order  to  facilitate  matching  the  boundary  conditions  for  the  Debye 

potentials  on  the  surface  r  =  a,  (120)  is  written  in  the  form  of  (116), 
realizing  that  the  function  Xn(^2r^  ^ecomes  infinite  at  the  origin 
where  the  field  (120)  does  not: 

P„(l>Os  «)  £« <t  (121) 

°22> 
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;  , «  - 1 .  jEMrf  -  » 


r~J2s±-< 


The  unperturbed  incident  wave  is  now  well  defined.  3ecause  the  incident 
wave  propagates  along  a  coordinate  axis  (z),  parameter  m  =  1.  As  a 
result,  the  Debye  potentials  all  have  the  associated  legendre  poly¬ 
nomials  of  degree  1.  The  internal  particle  fields  must  be  of  the  same 
form,  without  Xn(k,r),  and  must  also  be  finite  at  the  origin.  Also,  the 
arbitrary  coefficients  must  exist  in  order  to  properly  map  the  conditions 
at  the  boundary; 

%  i"''  U,r>  •)  e«  (123) 

r  7T/  r  i”'1  dntitK,r)P.W(  toss')  SIN  4  (124) 

The  scattered  wave  must  vanish  at  infinity  without  entailing  an  infinite 
energy  and  the  function  ^n(kr)  contributes  just  this  property  while 
maintaining  the  proper  fora  for  IT (r).  Arbitrary  coefficients  are  a- 
gain  necessary: 

rO  =  '7*  4  1  n~'  (kzr)P„0)(coia)  us  0  (125) 

r1T*a-j[*  %Jn" 

The  boundary  conditions  at  the  interface  ensure  the  continuity  of  the 
tangential  coiroonents  of  S  and  H.  In  soherical  coordinates,  this  trans¬ 


lates  into: 


Eq  (r=a)  =  (r:a) 

Etf  (rza.)  =•  E  of  C  r;  a.) 
He  (r-.o.)  -  H©5  (r-a) 


027) 

028) 

029) 


H$  (r* a)  =  Hf(r=a)  030) 

where  ( 1 )  refers  to  the  sphere  interior  and  ( 2)  the  exterior.  The  mix¬ 
ture  of  terms  in  JTj  and  TTj,  at  f-.rst  glance  makes  it  difficult  to  apply 
the  boundary  conditions.  However,  by  forming  the  appropriate  linear 
combinations  of  (81)  -  (86),  the  boundary  conditions  can  be  recast  into 


an  equivalent  form  for  the  Debye  potentials.  Define: 

•*  =  Ja  (s«ie£a)  +  T? 

(131) 

3  =  7w(SINSHg)  +  -fjs  H# 

(132) 

o  =  -^H9  (sins  HuO 

(133) 

=  =  &£e 

(134) 

Then  the  boundary  conditions  (127)  -  (130)  become: 

ACl^(rseO  -  <n'(r  =  a) 

(135) 

BCl,(rsa.)  =  BW(l-s<L'i 

(136) 

C  0}Crs<O= 

(137) 

(r--o.)  -  i>^(  r-.o.-) 

(138) 

Note  that:  —  IT  ** 

(139) 

ir(i,=  irstr* 

(140) 

The  quantities  A  and  C  remove  rf^and  3  and  D  remove  rfl,. 

The  decoupled 

equations  then  result  in  boundary  conditions  on  the  Debye  potentials.  i 

Substituting  (SI)  -  (36)  into  (131 )  -  (134)  yields: 

rA*  041) 

rg  a  C&S  *■  Slfle  tyj  042) 

rc  =  cm  [K  rlTj,^a  £^[<^3  (143)  i 

rD  ,  Cos  O&fcrlQ *  S/M  0  ifa/X  r HQ  ^  r ITJ  044) 


For  equation  (135)  -  (138)  to  be  valid  for  all  0  ,  0  on  the  coordinate 
surface,  r  =  a,  requires  the  various  bracketed  arguments  in  (141)  - 
f 1 AA)  to  be  enuivalent  betueen  regions  f 1 )  and  (2).  Thus: 


Now,  using  (139) 


JL 

dr 


£r<J 


.0)  r 

if,  rtr. 


and  (140): 

(149) 

(150) 

=  tf|('r(tr‘*  ff,S) 

(151) 

=  Yff  r(tr;.ir;) 

(152) 

Equations  (149)  -  (152)  are  the  required  boundary  conditions  on  the 
Debye  potentials.  Substitute  (121),  ( 1 23 ) ,  and  (125)  into  (151): 


Since  all  terns  in  the  series  expansion  are  linearly  independent  of  each 

other  for  each  value  of  n  and  far  all  & ,  gl;  then  the  corresponding  terns 

in  the  series  must  be  eaual: 

'yC *  ■>. 


0) 


^i-cnrn(n,o.)=K  jr„M-  3" ^*9] 
C„  :  k£>c^  |9n^la)"a-i 


Ka 

The  quantity-— — ^  can  be  reduced  by  using  equations  (S7),  (39),  and, 
recalling  that  the  frequency  of  an  electronagnetic  wave  is  unchanged  in 


traversing  the  interface: 

kJ  yc! 

Applying  similar  arguments  to  the  remaining  boundary  condition  yields  a 

Kj  rAjlCg  _ 


arsing  the  interface: 

d  r,w  -x?  K?x“’  x?  <w  . 

■I  Y?  =  -  vTf  x «  K?  '  =  X?  =  <  w  =  1 


term  of  the  form: 


K, 


=  m 


>/m2  = 


m 


The  four  boundary  conditions  then  become: 


CiYnC*.®-)1  mfVj/Otj&l-a/i  (***)]  053) 

4yn'(«,“)=  nr,Crn'(K2a)-l,r,<]n(K!,a.)]  054) 

C„yn(K,a)=  £r„(K.«)-ao^fK.4I  055) 

an¥„  (K.o.)  =  C'tt(Kia-yhn  (,56) 

ITote:  Y‘(Kr)s^YCKr) 

This  set  of  equations  uniquely  determines  the  coefficients  afl,  bn,  cr) 
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and  d  .  However,  only  those  corresponding  to  the  scattered  wave  are 


important  here.  Define  the  following  parameters: 
04  =  k,CL  -  Z  =  277*  miQ'- 

^  ?*o 

/3  =  tf,o.  =  ^ 

A  o 

where  A*  is  the  vacuum  wavelength. 


(157) 

058) 


The  quantity  «x  is  called  the  dimensionless  size  parameter.  Solve 
for  cQ  from  (155)  and  use  the  result  in  (153).  Also  solve  for  dn  from 

(156)  and  use  the  value  in  (154).  lie  obtain: 

%  (« )  VS  (a)  -m%(4)  YJ,  M  _  ,  . 

an  =  fa  (« )  Yn'CA) -/»  r*  f/3)  faJC<* j 

b  s  )  -Vn  (d)  #•>  '(•*■)  ,  (160) 

a  7*5 - C/3J 3/  (<*-) 

All  the  field  vectors  are  now  uniquely  described  given  the  known 

parameters  m  and**.  This  data  enables  the  coefficients  a  ,  b  ,  c  ,  and 

n  n  n 

d^  to  be  computed  and  their  potentials  determined,  finally,  the  vectors 
are  obtained  from  (Si)  -  (36).  The  formal  solution  is  complete. 
Practical  light  scattering  observations  are  performed  in  the  far  field, 
or  wave  zone.  This  condition  is  such  that  k2r  >J  n  where  n  is  the 
Ricatti-Pessel  function's  order  parameter.  In  this  approximation,  the 
previous  relations  become  simpler.  The  asymptotic  expansion  of  the 


Hanker  function  becomes  CO  • 

CntO  e-t«tr 


(lol ) 
(162) 


The  fields  become  transverse  due  to  the  radial  component  decay  as 
(A/r)  compared  to  the  (X/r)  dependence  for  the  ncn-radial  components, 
finally,  using  (Si)  -  (56)  with  (125),  (126)  and  (lol),  ( 1 62)  for  the 
scattered  wave  only  (  the  % (k^r)  tends  to  zero  for  large  h^r  with  re¬ 
spect  to  5fl(::.r)  as  does  V*/ )  yields: 


caused  by  an  incident  nave  of  intensity  cos  pfgolarized  parallel  to 
the  scattering  plane. 

The  ohase  difference  that  exists  between  these  two  oarts  cf  the 


, ’fl’/C  *!  P* 
*ijn  • . o. »  v  < 


,  r  Re  (Sj  Xm  )  ~  R&  lm  ^  i) 
tan^  =  fte(S1)Re(S«)+Im(S,)I«rSa) 


(173) 


Unrolarized  incident;  light  gives  rise  tc  a  scattered  bean  of  intensity: 

x*  r-  ■  \ 

074) 


Ju  =8i rV  (L*+i0 


with  degree  of  polarization  given  by: 


P  = 


<-|  ~  <•  2 
L  i  +  1 1 


( i  vO 
^ 1  •  ^  / 


Fixing  one 


.e  scattering  plane  m  space  to  ce  tne  yz-p_ane,  aii  co- 


servatior.s  by  definition  take  place  in  the  scattering  plane.  The  in¬ 
cident  electric  vector  has  its  direction  in  the  xy-plane  at  an  angle  X 
tc  the  y-azis.  The  x-axis  is  designated  vertical  and  the  yz-plane 
horizontal.  This  is  shorn  in  figure  3.  Thus,  the  vertical  (V)  and 
horizontal  (•{)  components  of  the  intensity  in  the  scattered  wave  for 
unit  incidence  become: 

<L 


^  ilStn%X 

I©  -  £  »/*)  r  lZ  C0S  *  X 

The  Stokes  parameters  are: 

S©  =  i/V*rl  fc  S/a/  ZX  V-  /'2  cos2  x) 
S,  =^T  g/  sw2X  -  /* 

^2.  S/A^  .X  X  C05 


076) 
077) 

(172) 
(179) 
(130) 

> j  -  fiirV*  v  Lt 1  z  **"’  (I'-i) 

’•Jhen  discussing  cross  sections,  one  is  speaking  of  the  ancunt  of 
area  that  the  particle  effectively  presents  tc  the  incident  bean  by  its 
ability  to  spread  out  the  incident  energy.  The  conplex  index  of 


| 
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1 ;  rwnpff 


refraction  includes  absorption  as  well  as  scattering.  Define  the  ex¬ 
tinction  cross  section  as  the  sun  of  the  cross  sections  for  scattering 
and  absorption  : 

(182) 

By  integrating  the  real  part  of  the  tine  averaged  Poynting  vector  for 
the  total  ext  .nal  field: 

033) 

the  total  outward  flow  of  energy  is  obtained.  ?or  unit  incident  energy, 
the  various  energy  losses  represented  by  the  integral  resulting  fron 


above  are  the  desired  cross  sections: 

(184) 

ni  (2n+,)  l  Re  + bn)  J  (185) 

The  efficiency  factor  is  defined  as  the  cross  section  per  unit 
2 

geometrical  area  (ffa  )  of  the  particle: 

TFa*  (186) 

Terms  of  the  formal  X  become sincere  =  2TTo./x  •  Thus: 

Qsco.  =  £  (Zn*i)  f  I O-rrl %/ 4>ni2?  (187) 

Qext  =  C^n+0  jRe  C&n+b«V}  (iss) 


Note  that  ^sca.  and^gjtfc  are  independent  of  the  state  of  polarisation  of 
the  incident  wave. 

Kerker  Q2J  et  al.,  gives  an  excellent  review  of  the  difficulties 

involved  in  the  numerical  computation  of  the  previously  listed  scattering 

functions.  The  functions  needed  to  be  evaluated  are  a, 

n7  n*  ”n'  Jn 

Some  discussion  of  the  equations  is  given  in  the  Appendix. 

Using  approximations  for  small  arguments,  the  fourth  power  law 
attributed  to  Rayleigh  scattering  can  be  directly  obtained.  The  first 
three  scattering  coefficients  are: 
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a, 

f.mlzL')  oc-1 
(.  +  z  J 

(169) 

0-2  S1s  i 

U/»t+3; 

(190) 

k|  “  WS  L 

(*>}*-/) 

(19 1) 

plus  terms  in 

higher  powers  of  «*% 

For  the  Rayleigh  criterion****!. 

only  contributes.  From  the  Appendix: 


IT,  (cos  Q')  z  l 
T,  (cos  9)  z  cos  e 
(165),  (166):  a.TT,  (cos  ©) 


(170): 


Sife)1  |  a.,  T,  (cos  e) 

;,co  =|  i«,r 


-J7  I  I X  CoS2  G 


For  unpolarized  incident  light,  for  example,  sunlight; 
(174):  I,  =^r 
(139):  Xu  ^ 

(157):  =^  (2-7rr)i’|P^r(M-c°s'le:i 


=>Xu  =  ff^  j^P(ltC05J©) 


(192) 

(193) 


(194) 


This  equation  accurately  predicts  the  blue  sky  phenomenon  as  the  shorter 
wavelengths  scatter  much  more  than  the  other  visible  wavelengths  in 
daylight.  Equation  (175)  yields  the  observed  degree  of  polarization: 


Id, 


a,lT-  la.pcps*© 
CL,!4  +  l<L,ta  Cos'  G 


?■>,  1 1  cos  9  _  (195) 

l  -f-cos  x  © 

A  plot  of  this  function  on  a  polar  graph  indicates  a  maximum  of  1 , 
complete  polarization,  at  a  scattering  angle  of  90°.  This  has  long  been 
observed  in  the  day  sky.  The  error  associated  with  the  truncation  after 
n  =  1  is  5%  if  a/x  is  kept  0.05.  This  limit  is  the  required  Rayleigh 


criterion, 


v/here  I0is  the  incident  intensity  to  the  volume  at  x-0.  This  last  e- 
quation  is  known  as  Bouguer' s  law*  [23}  It  describes  the  degradation  of 
the  incident  bean  as  it  traverses  a  material  path.  In  order  to  express 
the  scattered  intensity  into  a  solid  angle  dfL,ue  recall  that  the  in¬ 
tensity  scattered  per  particle  into  solid  angle  dXl  was  given  by  equation 
(200).  Let  the  parameter  r  in  that  equation  be  replaced  by  R  as  r  now 
describes  the  particle  radius. 

(215) 

As  before,  the  number  of  particles  in  the  cell,  II,  applied  to  equation 
(215)  yields  the  total  scattered  intensity  per  unit  radii  per  volume; 


then,  for  all  particles  of  radius  r  between  r^  and 

=>  X  =  X  (  m,  ©,  $0 

For  unpolarized  radiation: 

i  ^  (Li  +  Lz) 


(216) 


(217) 


where  i^  and  i^  are  given  by  equations  ( 1 65)  -  ( 1 63)  and  ( 1 70 ) .  For  the 
polarized  radiation: 

J£(c)s,nV+  u  cos  V)  (21s) 

Equations  (212)  and  (216)  are  the  specific  results  of  the  theory 


used  to  investigate  aerosol  parameters  and  will  be  discussed  in  the  next 


two  sections  of  this  paper.  First,  some  comments  on  the  collective 
scattering  phenomenon  are  in  order.  It  is  well  known  that  the  more 
ordered  the  scattering  centers  become,  the  more  coherent  the  scattering 
due  to  the  phase  relations  between  sites.  P  This  leads  to  a  reduced 
backscatter  i/ith  corresponding  increase  in  forward  scattering.  In 


crystalline  substances,  there  is  essentially  no  backscatter,  when  off- 
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. . 


resonance;  and  the  incident  bean  merely  undergoes  propagation  velocity 
changes.  The  condition  met  by  a  gas  at  atmospheric  pressure  or  less, 
where  the  number  density  is  low  on  a  radiation  wavelength  scale,  ensures 
that  the  scattering  is  incoherent,  allowing  the  summing  of  intensities 
from  all  the  scattering  centers.  The  Mie  theory  for  a  collection  of 
particles  requires  a  random  distribution  of  scattering  sites  to  achieve 
this  incoherent  intensity  accumulation.  Quantitatively,  the  near- field 
interactions  between  the  particle  may  be  in gored  if  the  average  distance 
between  sites  is  at  least  ten  times  their  radius.  This  criterion  is 
satisfied  even  for  very  dense  fogs  and  the  independent  scattering  model 
computations  are  valid  for  practical  tropospheric  scattering.  However, 
significant  number  densities  wil1  lead  to  a  large  flux  of  previously 
scattered  light  in  the  region  and  add  multiple  scattering  effects  to  the 
received  signal.  A  quick  check  on  the  validity  of  the  single  scatterer 
model  is  to  verify  linearity  between  the  observed  effects,  I  and/3,  and 
the  number  density  of  scattering  particles.  If  the  relation  is  non¬ 
linear,  one  must  suspect  the  existence  of  coherent  and/or  multiple 
scattering.  Alternatively,  if  the  relationship  is  linear,  it  is  certain 
that  the  independent,  single  scattering  conditions  are  present,  A 
measure  of  this  condition  may  be  obtained  from  equation  (214).  {24]  Lei 

Tx  s  -jf  (3  (s)  ds  (219) 


be  the  optical  depth  of  the  medium,  then: 

-rx 

I C*)  = 

The  conditions  are: 


(220) 


(transmission  >S0%)  -  single  scattering  predominates 
0.1$  r*<  0.3  (transmission  75  to  95%)  -  multiple  scattering 
effects  may  be  removed  via  correction  to  the 


IpS  pS 
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single  scatter  model 

Tx  (transmission  <75%)  -  the  complete,  non-ana- 
lytic,  approximate  multiple  scattering  theory 
applies 

Virtually  all  aerosols  in  the  troposphere  satisfy  the  first  criterion. 

The  literature  abounds  with  material  containing  graphs  of  the  Kie 
scattering  functions  for  many  specific  cases.  The  Appendix  outlines  one 
such  case  for  a  specified  set  of  condii  13.  However,  some  general 
information  can  be  reduced  from  this  plexor  a  of  data.  The  amount  of 
forward  scattering  is  about  equal  to  that  of  the  backs cattering  at  a 
size  parameter  (o*)  equal  to  0.01.  Asciis  increased,  the  amount  of  for¬ 
ward  scattering  also  increases  with  a  like  decrease  in  the  backscatter* 
At«<  -  1,  the  forward  scatter  is  on  the  order  of  100  times  that  of  the 
hackscatter.  Increasing®*  still  further,  there  exists  a  unique  scattering 
maximum  with  the  existence  of  smaller  subsequent  maxima  and  minima. 

The  scattering  function  increases  smoothly  as (Rayleigh)  up  until 
a/0.5.  Thereafter,  the  structure  is  highly  dependent  on  the  complex 
index  of  refraction.  The  large  oscillations  in  the  behavior  as**?  is 
varied  become  smoothed  as  the  absorption  index  increases.  As  the  num¬ 
ber  density,  as  a  function  of  the  particle  radius,  takes  on  more  of  a 
tailed  structure,  there  is  a  further  smoothing  in  the  scattering 
functions  as  small  differences  tend  to  be  averaged  out.  Conceptually, 
one  has  at  hand  either  an  actual  particle  distribution  or  a  specified 
model  of  such. 

Given  il(r),  for  each  r,  we  can  compute  the  required  (*•,©)  and  sum 
the  total  intensities  using  H(r)  as  a  weighting  distribution.  This 
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total  1(^,0)  then  appears  as  the  smoothed  function. 

lihe  the  scattering  function  publications,  there  are  also  a  large 
quantity  of  analytical  number  distribution  models  in  use.  The  verifi¬ 
cation  and  application  of  these  models  has  been  the  primer;/'  purpose  of 
much  of  the  aerosol  e:cperiments  conducted  in  the  past  20  years,  and  the 
nor!:  is  continuing  today.  As  computers  become  faster,  real  time  inver¬ 
sion  analysis  ’./ill  lead  to  a  more  accurate  determination  of  the  "true” 
number  distributions.  Two  such  models  shall  be  mentioned  here.  Some 
natural  aerosols  may  be  described  \-ith  a  tue-parameter  function  utilising 
a  multiplicative  scaling  coefficient,  c,  and  a  shaping  parameter /V,  as 
presented  by  Jungs  03  :  _  ,  . 


ol  r 


(221) 


d:?(r )  , 


’./here  is  as  previously  defined  in  (210)  and  V  may  vary  from  about 

a  r 

2.2  to  4»C«  As  an  example,  for  silicates  in  the  O.l^n  to  10.0  un  range, 


one  may  use: 

T,  =  O.04  yum  ?  rz  -  I O.o/t/n  ^  V  =3  . 

Another  such  model  is  in  reality  a  set  of  models  that  includes  the  Jur.ge 
model  as  a  special  case.  Deirmendjian  [iQ  has  proposed: 

Hr)  =  Mkl  „  ar"‘e-brT  (222) 


(222) 


(note:  the  used  here  is  different  from  the  size  parameter.)  The 
versatility  of  this  approach  is  the  presence  of  four  adjustable  parameters 
(a,©**,  b,  y)  allowing  for  a  wide  range  of  experimental  data  to  be  fitted. 
Three  specific  cases  of  this  model  are: 

Haze  C  (continental):  -P  ~  f  ^ 

-b'Tr" 

Haze  11  (maritime)  :  r  ~  re 

L  -1.5  r 

Cloud  : 
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Host  of  the  other  proposed  model  distributions  utilise  some  form  of 
equation  (222)  as  their  basis  and  canned  routines  for  their  computer 
modelling  are  readily  available. 

Two  other  effects  that  have  been  omitted  in  the  present  theory  are 
Doppler  shifts  and  Brillouin  scattering.  The  presence  of  random 
(thermal)  acoustic  naves  in  an  atmospheric  cell  ’./ill  add  a  degree  of 
order  "hereby  Brillouin  scattering  of  the  incident  nave  can  occur  off 
of  the  resultant  density  fluctuations.  Bo ’./ever,  this  effect  is  extremely 
■./cal:  for  gaseous  materials  and  the  degree  of  coherence  added  to  the 
structure  of  the  gas  is  generally  undetectable.  Similarly,  this  motion 
of  the  gas  will  cause  Doppler  broadening  of  the  wavelength  of  the  in¬ 
cident  radiation  adding  to  the  natural  line  uidtn  of  the  source.  Iven 
for  the  sharpest  sources,  lasers,  the  natural  linev/idths  are  on  the 
order  of  several  MBs.  '.Jhcroas,  for  almost  all  atmospheric  conditions, 
the  expected  Doppler  shifts  are  on  the  order  of  a  feu  tenths  of  a  Illiz. 
Therefore,  the  shift  is  generally  lost  in  the  natural  iineuidth  and  the 
relative  notion  between  the  detector  and  the  atmospheric  cell  may  be 
neglected. 


III.  LIDAR 


The  sophistication  of  the  laser  systems  available  today,  coupled 
with  the  fast  computer  processing  that  can  cheaply  be  obtained,  has 
singled  out  the  laser  radar,  LIDAR,  system  for  increased  use  as  a  tropo¬ 
spheric  experimental  device.  The  LIDAR  should  have  the  ability  to  per¬ 
form  a  wavelength  scan  with  selectable  bandwidth  receivers.  It  is  the 
atmospheric  attenuation  of  the  laser  pulse  and  the  introduction  of 
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random  phase  and  amplitude  changes  that  combine  to  form  the  scattered 
pulse,  or  echo.  The  subsequent  inversion  of  this  data  then  allows  the 
experimenter  to  assess  the  optical  characteristics  of  the  atmosphere 
on  a  real-time  basis,  with  the  system  geared  to  perform  azimuthal  and 
elevation  scanning  as  well  as  range  gating,  one  can  then  map  out  the 
atmospheric  optical  behaviour  within  the  visible  hemisphere.  The  two 
most  important  atmospheric  parameters  to  be  considered  for  a  pulsed 
LIDAR  system  are  the  extinction  and  scattering  coefficients.  The  total 


extinction  coefficient  may  be  ’,/ritten: 

@  *  /3  A  t  /30  (223) 

where  ft  -  Rayleigh  (molecular)  scattering 
H  -  Mie  (aerosol  +  other)  scattering 
A  -  Molecular  and  Aerosol  absorption 
0  -  Other  processes  (Raman,  fluorescence,  etc) 

The  other  processes  taken  together  are  typically  two  or  more  orders  of 
magnitude  weaker  than  the  Rayleigh  effect  and  tend  to  get  lost  in  the 
signal  noise. 

Rayleigh  scattering  is  obtained  from  the  Mie  solution  assuming  oC«l ; 
using  equations  (134)  ana  (1S9)  -  (191): 

lrR(sco-')  *hr 3  (V)  (224) 

Rayleigh  scattering  for  air  molecules  assumes  m  is  real  and  m  #1 ; 

*  '  *#■«/*> 

(r^sca.)  * -a  (m»- 1)1  V 1  (225) 

A. 

q 

where  V  =  (4/3)fTa  ,  the  volume  of  one  molecule.  Let  N  -  l/V,  the 


(224) 


molecular  number  density,  then: 

SF 


(226) 
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where  u  =  (&¥3§)/ (&-7<f) ,  the  depolarization  factor,  required  (227) 
to  account  for  the  anisotropy  of  air  molecules,  for  which  -  0.035*  [j&. 


The  quantity 


g(©,$0  a  O.0b07i  (1+0.932  oos2©) 


(228) 


also  accounts  for  the  anisotropy.  Let  (°Cu  be  the  air  molecule  number 


density  (vice  21 ),  then: 


/3«  =  1°*.°* 


(229) 


’where  it  is  assumed  all  air  molecules  have  approximately  the  same  radius 
of  10"4,i«n. 

The  Mie  extinction  coefficient  is  dependent  on  particle  radius 
through  equation  (212)  whose  only  unknown  is  dli(r)/dr.  The  unique  sharp 
wavelength  feature  of  the  LIDAR  allows  one  to  work  well  away  from  the 
characteristic  absorption  lines  of  the  sample.  Conversely,  it  also 
allows  one  to  type  the  scatterer  composition.  For  example: 

Laser  1:  \  on  an  absorption  (resonance)  line  of  a  particular 


species; 

=>  /3.  =  r  (3fl,  (230) 

where  (3^  is  the  resonant  absorption  extinction  coefficient  wiiich  is 
strongly  dependent  on  X  with  tabulated and 

K  =*Af>'  (23,) 

The  (°  is  the  number  density  of  the  species. 

Laser  2:  just  off  the  resonance; 

& 

@z  ~  @RZ  +  +  @42  ^3^ 

Due  to  the  relatively  weak  X  dependence  for  wet 

(233) 

For  an  analysis  cell  a  distance  R  away: 


60 


(2U)  I  (R)  =  I.  e'C*Cs)is 


and  back  at  the  receiver:  R 

l'(0),X(R)e  -f^CsUs 

X,'  X0,  e’%R/3A(s)ds 

Then:  -~rr  -  =—  e 

-i-2  -u<>2 

Thus,  the  t®1  for  the  particular  species  nay,  in  principle,  be  found  fron 


(234) 

(235) 


this  equation  using  the  inversion  techniques  outlined  in  the  next  section. 
For  the  typical  nagnitudes  of  the  considered  cross  sections,  see  Table  II. 

TABLE  II 


Process 

0"  ( cm  /  ster) 

Hie 

10“ 27  to 

]G~° 

Rayleigh 

10“27 

Raman 

10“ 30  to 

io-29 

Resonance  Raman 

10"  23 

Fluorescence 

£10“16 

Similar  to  radar  in  operation,  the  LIDAR  has  a  performance  equation 
’./hereby  the  system  characteristics  and  the  appropriate  atmospheric  para¬ 
meters  vi'Cih  contribute  to  the  signal  are  related. [23)  Consider  a 
pulsed  LIDAR  operating  in  nones tatic nw.de  (receiver/ transmitter  co-located) 
being  used  to  analyze  a  specific  atmospheric  cell  via  its  optical  prop¬ 
erties.  The  problem  geometry  is  outlined  in  figure  9.  In  this  figure: 
hx  -  cell  height 
R  -  slant  range  to  cell 
a  -receiver  solid  angle 
-  transmitter  solid  angle 
L  -  cell  length  "gated”  via  receiver  on-time 
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Monostatic  laser  racar 


laser 

transmitter 


laser 

receiver 


4 


(exaggerated) 


A  -  cell  area  seen  by  receiver 
Ax  -  cell  area  seen  by  transmitter 
The  following  derivation  shall  consist  of  counting  photons  at  each 
stage  along  the  bean  path;  converting  to  intensities  has  been  previously 
discussed. 

1.  II.  -  number  of  photons  transmitted  cer  oulse 

O  A 

\JL  -  photon  frequency  (transmitted) 
hVj.  -  energy  per  photon 
H.h^  -  energy  per  pulse 

"C  t/  w 

Pj.  -  laser  output  power 
£\  -  pulse  width 


-OL 

h  Vt 


(236) 


2.  Per  nonostatic  LID  AT.,  the  scattering  angle  is  ICO  degrees  (Tf  } 
for  backscatter.  Assume  that  all  hacks cattered  photons  are  received: 

=,>  nR  ^  rit  (237) 

The  cell  volume  is  determined  by  the  IZDAR  timing  sequence  and  is: 

V  =LP* 


?ne  solia  angle: 


n  *  fl/R! 


V  =  LAtR 


•v  j 


(239) 


(240) 


The  pulse  width  constrains  the  maximum  system  resoiu-cior.  by  limiting  the 


minimum  ceil  length  achievable: 


Lm,„  r  y2(crt) 


(241) 


where  c  is  the  speed  of  light  in  the  medium. 

This  condition  assumes  that  the  pulse  tail  reflecting  off  the  front  face 
of  the  cell  does  not  interfere  vith  the  pulse  front  reflecting  off  the 
rear  of  the  ceil. 


W  ,„  ,  «-  ««m  .  ,»,  »,f  t.  '-  '”  !  "  "  -i 
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3.  -•-  -  number  of  photons  incident  on  the  cell 

(3^  -  atmospheric  extinction  coefficient  for  the  given  ^(V*} 


transmitted 

hlL  sNte. 


(242) 


however,  to  account  for  bean  forming  requirements,  assume  an  optical 


system  efficiency  of  ‘y|i: 

=>  NL=ntNte'^fs)ds 


(243) 


4.  The  backs  cat  terir.g  coefficient  nay  be  defined  in  accordance 


with  equation  (212): 

(°  a  (ft)-  (3  (ft)  - J\ (r»,  \n  r,  6  =  IT,  0  <*r  (244) 

where  ^crCu)  is  the  usual  notation.  This  equation  represents  the  back- 
scattered  intensity  per  unit  incident  intensity  per  unit  solid  angle  per 


unit  path  length.  The  signal  returned  Iron  the  cear  is  due  to  this 
process  alone.  3cuguersT  law  (214)  states,  for  (°<3(ft)  assumed  constant 


throughout  cell  length: 


X-Ioe 


-(Oor(ff)L 


(245) 


where  I  is  that  uorticn  of  the  incident  beam  that  masses  through  L  un¬ 


affected  .  The  scattered  portion  is  I  =  I„-  I  : 


f/s  =  NL(l-e-^Wl) 


(246) 


where  N  -  the  number  of  mho  tons  bad-:  scattered  ar.d  a  ratio  of  intensities 
s 

is  the  sane  as  the  ratio  of  the  number  of  photons.  Typical  tropospheric 
aerosols  yield: 


(O  era OL  «  I 

(i_  (°cr(n)i  -  0.02,  overall  error  IT.) 


(247) 


The  restriction  this  mlaces  on  1  is  a*,  upper  limit  of  several  hundred 


r“}^T  • 


/OO-ML 

=>  Ng  =  tiL/°tr(rr)L 

The  urine  or.  ”  indicates  the  number  cf  oho  tons  cer  solid  angle  t 
s  • 

been  backs cattered.  Ocnoining  equations  (243)  and  (24?): 

N^s  =  tJt  %/°<r(r)L  e%P  f-jTfy  W  ds  } 

5.  '"or  the  echo  must  pass  through  the  sane  propagation  path 
return  to  the  receiver  as  it  did  cn  its  way  to  the  ceil,  however 
extinction  coefficient  nay  be  different,  (3  }  due  to  the  accunulat 
“'recesses  like  kansn  scattering,  Irillcuir.  scattering,  loonier  shnrt. 


9  j  oi. 


-  nunoer  o"  riiooons  at  receiver  rirror 

Nm  =  ^nre  ~-{ferCs^s 

(239),  (250):  n  ,  1 

=  ^7 t/°o-Crr)LPrexp  f~/a  fate+'3r(s}jdjl 


(251) 


(252) 


6.  The  optical  systen  ’.dll  collect  the  echo  photons  received  at  an 
efficiency  of  ^  _  and  convert  then  to  an  electrical  signal  via  a  detector 
t.’ith  a  quantum  efficiency  of  7  : 

^f=  NrwV?rV?<j,  (253) 

where  N  -  number  of  ohetens  received  and  counted  oer  transmitted  oulse, 
r  7 

due  to  backscatter  fron  the  atmospheric  ceil  of  interest. 

Thus:  /•  /-fly-  ^ 

Mr=  ^%Par(rr)Lfi^r^^PLJ0  dSj  .  . 

ft2 

This  equation  nay  be  recast  in  the  fora  of  intensities  if  the  photon 

population  returned  is  large  enough  that  individual  photons  need  not  be 

counted.  Let  I)  be  the  ootical  systen  aperature  (Dx  -  transmitter;  D  - 

»  r 

receiver),  then  the  intensity  becones: 


mv  ' 


X  =  NhwD 

Reciting  equation  (254):  r  n  1 

lr-.  -  — ^ 


V03j 


(256) 


Thus,  the  important  parameters  for  the  LUSAR  analysis  of  an  atmospheric 
cell  are  the  scatter er  number  densities  and  associated  cross  sections. 

For  a  given  system  design,  one  has  available  the  Quantities  Ii} *), , 
A^,  D^.,  Dy,  and  can  select  the  values  for  %)* ,  1,  and  R.  If  one 
assumes  a  slow  moving  target  cell,  then  the  pulse  round  trip  tine,  ty,  is 


tr  =  ZR/c 
•>  R.=  (c-A)^r 


(257) 

(25S) 


Thus,  one  determines  the  range,  R,  by  turning  the  receiver  optics  on  a 
time,  t^,  after  the  pulse  is  emitted.  In  order  to  avoid  the  problem  of 
range  ambiguity,  like  radar,  one  must  ensure  that  t^  be  less  than  the 


period  between  pulses: 

tr  <  \JPRF 

where  FRF  -  pulse  repetition  frequency 
The  unambiguous  range  limit  may  be  written; 

K  =  (c  lz)/PRF 


where  R  -  unambiguous  range 

u 


(259) 


(260) 


The  cell  length,  L,  is  determined  by  the  time  that  the  receiver  is  turned 
off  again  a  tine  At  after  ty.  For  the  entire  cell  to  be  scanned,  the 


mdse  must  traverse  L  twice  in  At: 


U  U/2)At 


and  the  receiver  must  be  turned  off  at  time: 


t  =  tr  i-  At 


(261) 


(262) 


One  can  then  select  L  and  R  by  programming  a  timing  sequence  for At  and 
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t  •  Finally,  the  interpulse  period,  T,  vliich  is  a  systen  dsadtime  be¬ 


tween  pulses  available  for  processing,  ceccnes: 


T-.  0/psf)-  rt 


(263) 

Now  one  measures  I  and  'l)r  and  applies  equation  (256)  to  detemine  the 
desired  quantities  (°cr  (ft)  and  (3 60.  This  inversion  is  covered  in  detail 
in  the  next  section.  Note  that  vhen  the  equation  for  ts(rn,y.)r,B,qf)  iron 
(202)  is  used  to  replace  the  jO  in  (244),  then  equation  (256) 

for  I  has  the  correct  fom  as  given  by  (21 6)  for  the  scattered 
intensities. 

Thus,  a  model  has  been  constructed  whereby  the  system  parameters  are 
related  to  the  cell  unknowns  in  a  definite  manner.  But  first,  some  con¬ 
sideration  must  be  given  to  systen  noise.  During  the  tine  the  receiver 
gate  is  open, At,  the  systen  will  collect  and  count  the  background  sky- 
photons  as  well  as  the  signal.  The  continuous  "dark  current"  in  the 
detector  yields  an  electrical  noise  as  a  result  of  dark  counts.  These 
counts  arise  from  photoelectrons  being  emitted  from  the  detector  photo- 
cathode  surface  even  without  an  illumination  present.  Each  of  these 
noise  sources  nay  be  reduced  by  proper  design  of  the  equipment  and  ex¬ 
periment.  In  the  case  of  a  detector,  its  dark  current.  L,.  is  either 
specified  or  nay  be  measured  in  the  laboratory: 

XD=NDeG  (264) 

where  -  number  of  dark  counts  per  sec 
e  -  electron  charge  per  count 
G  -  detector  count  gain 

The  number  of  dark  counts  detected  as  pare  of  the  return  signal  in  a 
tine  At  is  IT-A t,  ana  : 


Al  jv  "  /(«(?)  (265) 

Let  Ng  represent  the  background  count  rate  such  that  II^At  is  the  num¬ 
ber  of  background  counts  received  per  laser  pulse.  Assume  a  background 
spectral  brightness,  \l\,  that  accounts  for  all  the  natural  and  artificial 
background  sources.  This  represents  a  power  density  per  solid  angle 
in  the  background  sky.  Then: 

=  (*H^r)  ~  A  (266) 

where  ^  v^fyCLri Dr  -  previously  defined 

V  -  frequency  of  background  photons  (atX) 

A  X  -  receiver  banawith,  centered  about  X 
Let  1;^  be  the  total  number  of  noise  counts  per  laser  pulse;  using 
equation  (261)  for  At: 

Nn  =  (ZL/c)(Nd  +  Nb)  (267) 

Due  to  the  pulse  repitition  frequency,  PR?,  the  noise  count  rate: 

M  =  (PRFjNr,  <2®> 

where  II  -  noise  count  rate 

The  statistical  nature  of  noise  yields  an  effective  RMS  noise  count  rate 
of  Vn  .  The  signal  has  a  count  rate  of: 

Sr  =  (PRF)Nr  (269) 

where  -  count  rate  from  signal  echo  alone 
Thus,  the  total  counting  rate  is: 

5=  PRF (Nr  +N)  (270) 

Talcing  the  RMS  of  the  signal  plus  noise,  one  can  form  the  signal-to- 
noise  ratio: 

S IH  =  Sr/J5~  (271) 

Assuming  a  one  second  time  averaging: 


6S 


4.  A  sn?21  ilr  requires  a  lower  bean  divergence  at  the  transmitter 
and  an  optical  matching  between  the  receiver  and  the  echo  bean  pattern 
(another  laser  capability). 

5.  Without  any  noise,  there  still  exists  a  finite  (3/ll)n  which 
represents  the  "shot  noise"  statistical  nature  of  the  signal.  It  must  be 
emphasised  that  ail  the  quantities  in  equation  (276)  represent  the  mean 
of  the  indicated  variables. 

6.  if  Nr  «\^n  =$  S/N  ~  Nr  (277) 

~  VjSL  (273) 
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Thus,  the  addition  of  (3^  to  the  total  (3 r  is  on  the  order  of  a  feu  per¬ 
cent  and  any  quantitative  measurement  must  take  into  account  both  the 
aerosol  and  molecular  scattering.  The  turbidity  is  defined  as  the  rati; 
of  the  total  scattering  intensity  to  the  Rayleigh  scattering  intensity: 


T, 


(279) 


=  Xt/Xr 

Some  researchers  attempt  the  measurement  of  this  quantity  independent  of 
the  intensity,  then  use  the  result  to  deduce  the  aerosol  I,..  Another 
method  assumes  that  the  aerosol  concentration  above  a  certain  height 
(e.g, >  30  km)  is  sufficiently  low  to  allow  it  to  be  neglected.  Then 
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v  v.  A-; 


equation  ( 256} : 

XR  ~  t%  OrW 


(2:0) 


r* 

Text,  one  computes  Ip  from  (256)  assuming  the  conditions  of  the  U.  3. 
Standard  Atmosphere  and  equation  (226).  “This  is  the  expected  echo 
intensity  due  to  molecular  scattering  alone.  The  ratio  between  the  mea- 
sured  and  the  expected  Ip  yields  a  scale  factor,  b,  that  corrects  the 
standard  atmosphere  for  the  local  conditions.  Using  this  factor,  one 
assumes  a  density  distribution  for  the  local  atmosphere  the  sane  as  that 
for  the  standard  and  commutes  the  exoected  value  of  at  the  cell's 
location.  This  intensity  is  then  scaled  using  b.  The  measured  value  of 
I_  is  then  corrected  for  this  I~(b)  and  the  result  becomes  I„,  the  Mie 

JL  it  ti 


intensity: 


XM  =  XT  -  bXj 


(231) 


?or  other  means  of  separating  the  molecular  and  aerosol  intensities  from 
a  measurement,  the  reader  is  referred  to  the  work  by  Hall,  et.  al.  {jO 
This  Rayleigh- atmosphere  assumption  is  the  starting  point  for  the  in¬ 
version  process. 

The  Standard  Atmosphere  may  be  represented  via  the  follouing  equations 

temperature  T(°F)  =  59.31  ~  0.00357  C'F/ft)  b(Ft)  (2S2) 

pressure  1,022.  £xp(-4.IS7)Uo'50/fOh(tt))  (233) 

Mass  densi  ■fer  -  o.OO?.MZJei.p('3.3S^^‘60/Ft)K^})  (2S4) 

(slugs/ft*7) 

where  h  may  vary  from  0  to  37,000  feet.jjjj  The  more  precise  form  for 
these  relations  are  also  represented  as  canned  computer  routines. 

The  critical  tern  in  equation  (256)  is  the  arm(fl')  where  the  sub¬ 
script  M  refers  to  the  aerosol  distribution  alone.  The  extinction  com¬ 
ponent  in  the  exponent  of  e  in  this  equation  represents  the  scattering 
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and  absorption  cross  sections  integrated  over  the  entire  4lf  steradian 
sphere.  If  one  assumes  some  representative  model  for  the  aerosol  dis¬ 
tribution,  as  was  done  for  the  Rayleigh  case,  then  the  only  unknoim  in 
(256)  is  'The  two  averaging  processes,  integrating  aver  41T  and 

integrating  along  the  path  R,  reduces  the  sensitivity  of  the  exponential 
term  in  (256)  on  the  actual  aerosol  distribution.  After  completing  the 
inversion  process,  one  may  input  the  best  estimate  for  this  distribution 
into  the  exponential  term  in  (256)  and  then  repeat  the  inversion.  Ko:?- 
ever,  the  processing  time  involved  may  prove  to  be  the  final  constraint. 

The  LIDAR  experiment  yields  the  scattered  Mie  intensity  from  the 
aerosol  population  in  the  analyzed  atmospheric  ceil  as  a  function  of 
the  pre-selected  parameters  at  the  laser  (e.g.  wavelength  and  polarization 
scanning).  With  this  data  as  input,  the  system  is  prepared  to  perform 
the  inversion  routine.  Some  computational  needs  are  reduced  for  the  case 
of  laser  backscatter  as  the  scattering  angle  remains  fixed  at  ISO  degrees. 
This  backs cattering  condition  simplifies  several  calculations  as  out¬ 
lined  in  the  Appendix. 

IV.  MATHEMATICAL  IH VERSION 

The  assumption  made  throughout  this  paper  implies  that  the  experi¬ 
menter  need  only  apply  the  exact  Mie  theory  to  his  observations  and  per¬ 
form  a  simple  data  inversion  to  discern  the  aerosol  properties  of 
interest.  However,  this  inversion  is  by  no  means  simple.  Aerosol  prop¬ 
erties  may  be  separated  into  two  distinct  categories.  The  first  con¬ 
siders  those  parameters  that  are  explicit  in  the  previously  defined 
scattering  functions;  i.e,,  spherical  aerosols,  aerosol  sizes,  low 


X  -  incident  radiation  wavelength 
Q,  ft  -  scattering  angles 

The  large  number  of  unknowns  that  need  to  be  determined  will  ensure  a 
non-unique  inverted  solution  from  the  limited  sample  of  observa.tions_.ar- 
vailable.  There  is  also  the  added  difficulty  of  variance  in  the  sup¬ 
posedly  known  data  whose  uncertainties  at  best  generate  unbounded  noise 
in  the  final  results.  Therefore,  some  attempt  must  be  made  at  reducing 
the  initially  large  number  of  unknowns. 

The  dependence  of  the  intensity  on  the  size  parameter,  has  been 
discussed  previously.  Briefly,  the  non-monotonic  fluctuations  as  o<  is 


process  to  be  of  good  use  in  tropospheric  analysis.  These  marina  and 
minima  are  reduced  by  increasing  the  field  of  view  and/or  the  transmitt 
bandwidth.  It  is  to  LIDA?.'  s  advantage,  then,  mo  retain  this  complex 
structure  as  a  data  filter.  See  the  Appendix  for  a  sample  plot  of  ?{oC) 
The  complex  index  of  refraction,  n,  has  a  dispersion  relation  like 


The  less  than  1/  change  in  the  index  of  refraction  is  typical  for 
dielectric  media.  The  change  in  the  absorption, y,  is  even  less.  Stud 
have  indicated  that  the  disrersior.  in  n  is  balanced  by  the  randomness  o 


dispersion.  \23J 

The  angular  dependence  is  a  pre- specified  sot  of  conditions,  espec¬ 
ially  for  the  case  of  a  1T.DAR  backscatter  system.  This  system  consists 
of  a  linearly  polarised  wave  undergoing  only  bacl: scatter  observations. 
As  a  result  of  these  and  similar  considerations,  one  can  reduce 


+  ' 


their  radii.  By  the  use  of  pre-selecced  values  for  r,  one  only  has  to 
solve  for  the  distribution  functions  granting  that  the  radii  range  is 
properly  considered.  Therefore  the  parameter  of  urinary  interest  is  the 
scatterer  number  density. 


Define  the  Hie  function,  K: 

K  (  1Y\, 0)  =  ^277 'r)  i 

and  the  particle  size  distribution  function,  f(r): 


dr 


(257) 

(258) 

(289) 


The  scattering  intensity,  equation  (21 6),  becomes: 

IW  =  f  V(X,r)  ?(0 

-'r, 

where  x  represents  the  set  of  user  defined  parameters  (a,  X,  Q ,  0).  The 
value  I(x)  is  measured,  the  quantity  K(x,r)  is  computed  and  the  unkr.o-./n 
f(r)  must  be  solved  for.  Hois  relation  is  classified  as  a  Fredholm  in¬ 
tegral  equation  of  the  first  kind.  The  kernel  of  the  equation  is  K(x,r) 
and  I(x)  represents  the  transform  of  f(r)  utilizing  the  stated  kernel. 
There  exists  only  a  finite  set  of  measurements,  I(x),  over  which  the 

■i 

inversion  may  be  performed.  Let  the  subscript  i  represent  the  i' 
measurement  of  a  set  of  n  observations: 

(29-0) 


.tn 


leasurenexx^  ui  a  seu  vx  a  ooservav/ions: 

IW  =  3^0-  91  -  f  \  (*i,r)f(r)  dr 

’I 


The  notation,  g^,  has  been  introduced  and  will  be  used  later  ir.  the  paper. 
The  inversion  process  is  by  nature  an  unstable  one  and  this  may  easily  be 


demons  orated • 


Define: 


$F(xi£)=  pKttiiOE^*  C  smFr]  dr 
Jo 


(291 ) 


where  i  =  1,  2,  ...,  n  F-1,  2,  3,  ... 

and  C  is  an  arbitrary  constant.  Letting  r^  -  0  and  -  IT  in  equation 

(290): 
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rir 

s  j(xi)  +  C  lsin(Fr)K(*i,r)dr  (292) 

Since  the  sin(Fr)  acts  as  a  ‘weighting  function  for  the  kernel  on  the 
interval  (0,11),  as  F  increases  without  bound,  the  second  integral  tends 
to  0.  Thus: 


(293) 


One  can  aooroxinate  the  transfomed  values  of  f  with  the  function  S_  as 

*  *  rt 

» 


clcse  as  possible  by  talcing  F  large  enough.  Yet  the  jCsinFrj  function 
nay  be  made  arbitrarily  large  and 

f(r)  =  f(r)  +  CsinFr  (294; 

will  transform  as  f  (r)  with  little  to  no  resemblance  to  f  (r)  for  large 
values  of  C.  Therefore,  one  can  generate  transforms  that  come  as  close 
to  the  desired  values,  g(x),  within  any  desired  accuracy  as  one  nay  wist 
and  still  have  no  idea  what  the  accuracy  with  respect  to  f(r)  say  ce. 

The  technique  used  to  remove  this  instability  is  tc  require  the  ad¬ 
dition  of  a  constraining  relation  on  f(r)  other  than  equation  (290). 

This  nay  be  done  either  explicitly  or  imnlicitly.  The  proper  relation 
must  contain  the  desired  qualities  of  f(r)  requiring  a  high  degree  of 
a  priori  knowledge.  Out  of  the  near- infinite  set  of  possible  solutions 
for  f(r)  generated  by  the  inversion,  the  constraining  relation  will  act 
as  a  filter  to  obtain  the  unique,  or  near-unique,  solution. 

The  presence  of  measurement  error  in  the  g(x^),  larger  than  some 
quadrature  error,  allows  the  use  of  an  appropriate  quadrature  approx¬ 


imation  of  the  form: 


(295) 


for  k  intervals  of  the  argument  y  and  the  u.  are  the  quadrature 

J 

weighting  coefficients  with  increments  Ay-.  For  equation  (290): 
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9W=  \  W;:  KC^ijrj)  f  (ri)h  rj 

jit  J 

where  i  represents  a  measurement;  and  3  a  radii  interval. 


(296) 


Define;  for  n  measurements: 


1.  column  vector  G:  g.  =  g(x. )  i  =  1,  n; 

2.  coefficient  matrix  /I:  a..  a  w.  .K(x.,r.) 

ij  i*  3 

j  -  1,...,  1;  k  radii  intervals 

3.  column  vector  II:  II.  -  f(r,.)  Ar. 

3  d  3 

Then,  equation  (296)  nay  be  written  as  a  matrix  multiplication: 

G  =  AIT 


(297) 

(29S) 


(299) 


(3C0) 


The  value  IT.  represents  the  total  number  of  particles  per  unit  volume 


with  radii  between  r.  and  r .  +  A  r . ,  as  can  be  seen  via  (238): 

.  3  3  3 

Nj  -yrj  f(r)3r  s  f  (rp  [(rj  +  A  rj)-  rj]  (301) 

and  f(r. )  is  considered  constant  over  the  interval  j.  It  is  further 
3 

assumed  that  the  scattering  functions  have  only  small  variations  with¬ 
in  each  sub- interval  radii  range. 

The  problem  has  now  been  reduced  to  solving  for  the  k  unknowns  in 
IT  given  the  (k)x(n)  computed  values  in  A  and  the  n  measurements  in  G. 
It  must  be  emphasized  that  a  quadrature  approximation  will  introduce 
serious  error  in  f(r)  unless  the  quadrature  weighting  formula  is  of  an 
accuracy  higher  than  that  of  the  expected  measurement  error.  Let  nA 
be  the  number  of  independent  observations  made  such  that  nA  jC  n.  If 
k  =  rv.,  then  II  may  be  solved  for  exactly: 

(require  A  to  be  non-singular) 

=>  A"’A  =  f\  A'1  =  1  (302) 


where*  1  in  a  matrix  equation  represents  the  identity  matrix, 

(300):  AN  =  G  =>  N  =  A'!Gr ,  (K-nP 


(303) 


where  II  represents  the  interpolated  solution  vector.  However,  if  k 
K.  nQ,  then  ;I  nay  be  solved  by  minimizing  the  square  of  the  error  (least- 
squares)  : 

(require  A  A  be  non-singular  s>  A  :  transpose  of  A) 

(300):  (ArA)H=  AT£r  (304) 

this  step  forms  a  square  matrix  which  may  be  inverted, 

M  =  (ATA)  '  f\T£r  y  (K-fO  (305) 

That  this  last  equation  actually  minimizes  the  error  nay  be  visualized 
by  considering  an  appropriate  spatial  geometry  which  contains  the  vector 
G  and  the  vector  AN.  where  N.  ^  :!.  See  figure  10.  The  vector  Q— AIT^  is 
seen  to  be  representative  of  the  error  vector  between  G  and  AN ^ .  The 
minimum  error  occurs  when  this  error  vector  has  the  smallest  magnitude. 


tnus: 


(AN,)  -(0-flN,)  =o 


(306) 


(the  •  represents  the  vector  dot  product). 


Then  for  any  AN  in  the  range  of  A,  the  vector  AN.  is  still  perpendicular 


to  it: 


(AW.)-  (G-AN)  =  o 

Writing  the  dot  oroduct  as  a  matrix  nultioli cation  requires: 

(AW,)r(c,-flW)=  0 


(307) 


(308) 

(30?) 


=,>  !V(AtG-AtAN)  =  O  (30?) 

And  Iiy  r  0  since  it  represents  any.  II ^  vector  in  the  range,  thus: 

(»4)  ATflN=AT& 

In  the  work  that  follows,  when  writing  the  inverse  matrix  A-1,  if  the 
matrix  is  not  square,  the  quantity  (ATA)~1AT  is  implied. 

The  procedure  as  indicated  by  equations  (303)  and  (305)  are  in  large 


unsuccessful  as  they  do  not  contain  the  constraining  criterion.  3cne 
further  problems  concern  the  size  of  k. 

A 

1.  Small  k:  the  vector  I:  does  not  accurately  represent  II  due 

to  the  r.  set  sparsity  and/or  the  crudeness  of  the  quadrature. 

J 

A 

2.  Large  k:  the  components  of  I;  take  or.  ar.  oscillatory  be¬ 
haviour  irrespective  of  II  and  the  coefficients  in  A  become  particularly 
sensitive  to  round  off  error. 


-‘a*;  2^0^  g  -  y*l  T*,*  n  'n  00  0”  *  "•  ~  g  t  0"r  * 


*■?  a?  i  f  AT  5n*'  T  **'»'*  *" ^  1*7  * 

A—  -'•*  **•*“•  ws"*  V*;—  ••  %v-  |>J  • 

to ncJ(fl)  =  jflj  jR'1! 

/here  |.’j  —  ncm  of  a  matrix 


■  ~1C' 
\  ^ 1  ~  / 


ifll  -  max  s  I  i 

1  ■  i  l^j^K  i>i  |  CLi-j  I  :  maximum  ct’.umK  -  sum  (-*'  >  ; 

where  | a,..)  -  magnitude  of  the  component  a.... 

IflL  =  16^ n  I CLijl  ^  }  :  ~a::inuz  return  C312; 

(normally  used)  jAle  =  (^,  J",  ^  kiji  :  Luolizc-ar.  norz  (313> 

Lecause  A--  must  be  computed  and  will  contain  at  least  the  inaccuracies 
in  A,  the  conc(A)  will  be  costly  tc  compute  end  inaccurate  as  •••ell. 
nevever,  tne  test  is  tne  magnitude  oz  me  ccr.o'ii, .  n  large,  one  maun 
A  is  termed  poorly-conditioned  and  the  likelihood  of  a  successful  ir.vcr 
sion  decreases  rapidly  to  zero.  Corrections,  like  scaling  or  partial 
pivoting,  must  then  be  applied  to  A.  Ilote,  for  the  identity  matrix. 


cond(l)  =  1; 


cond  (A)  ^  ! 


U«  Li 


since  A  must  exist.  Other  orouerties  cf  the  no: 
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A !  ^  0 

l  At  =0  ‘iff  A  -  o 
|  /3  A  1  =  A  I A  I  y  ft-  >  0  a nd  cl  number 
Ifl+BI  ^  Ifll  +  i  Bl 
IflAl  <  1  Bl  I Al 


(316) 


(317) 


(31-3) 


(319) 


The  ccnd(A)  is  most,  useful  in  bounding  the  error  expected  in  the  computed 


solut^r'  and  relating  this  error  to  the  magnitude  of  the  residuals. 


v-  ■. 


\y  s  3 


(320) 


,:here  X  -  exact  solution  vector 


2  -  approximate  solution  ’rector 


3  -  vector  oi  xnoun  Quantities 


A  -  operator  matrix 


Tne  resiausus  are: 


R  =  3  -  AX 


(321) 


where  R  -  residue!,  vector 


3*2-  X 


(322) 


R  =  A3 


(  o  O  ^  N 
} 


Using  3  by  inverting  equation  (323)  and  the  norm  relations  (31 5)— (31?/ > 


one  nay  incorporate  the  cord  (A)  as  an  error  measure: 


JudM  IB!  4  ifl"  ^ 

A 


(324) 


Therefore,  the  relative  error  in  X,  contained  in  3,  can  be  as  great  as 


the  relative  residual  times  the  cond(A)  and  as  small  as  this  same  re¬ 


sidual  divided  by  cond(A).  Tor  large  cond(A),  the  residual  yields  little 


to  no  information  on  the  accuracy  in  X.  Equations  (320)-(323)  may  be 


used  to  define  or.  iterative  improvement  algorithm  (IIA). 


O' 


?>  %  h 


A  A 

1.  compute  X  by  X  -  A  3 

(321)  2.  compute  R  by  R  =  3,  -  AX 

A  A  _  1 

3.  compute  3  by  3  -  A  ?. 

4.  correct  X  bj  X  =  3  +  X 


(325:- 


(326) 

(327) 


If  the  system  matrix  A  has  errors  in  its  coefficients,  the  system 


may  then  become  ill-conditioned  even  if  the  residuals  are  small. 


A  =  A  + 


(32S) 


'./here  A  -  true  coefficients 


-  error  m  coefficients 


The  equation  actually  being  solved  is: 

AA 

a::  -  3 


(329) 


Then,  analogous  to  (324): 

-*,Vr  6  co',i>  (330) 

The  relative  error  in  the  residual  can  then  be  no  larger  than  the  rela¬ 
tive  error  in  the  coefficients  in  A  times  the  cond(A). 

Finally,  ar.  approximate  accuracy  check  may  be  obtained  from  Jj>3 : 

lE  or  HL  _p  . 

if  \%\  ^10  ,  (331) 

A 

then  X  is  probably  correct  to  p  digits. 

Equation  (300)  is  assumed  to  be  accurate  within  a  small  measurement 
error  and  equations  (303)  and  (305)  are  considered  insufficient  to 
remove  the  inversion  instability.  There  are  several  techniques  avail¬ 
able  to  remove  thses  instabilities  of  which  three  will  be  briefly 
covered. 

The  first  method  attempts  to  smooth  the  instabilities  by  one  of  tw; 
means.  Define: 

u(n)=(g--AN)t(g-AM)  +  v,  (N-N0)t(n-N;)  <72(0M)t(BN) 


(332) 


JHkW^Tv  CWWwWIIMWWWWWIJWI  '  PwaHre^wSS® 


or 


if  -  31T 


(352) 


in  matrix  forn  ’./here  3  is  a  (k)x(k)  matrix  defined  as: 

f  1-210 

0  1-2  1 

0  0  1-2 

0  0  0  1 

I 

3  = 


0  0  0 
0  0  0 
0  0  0 


1  -2  1 

0  1  -2 
0  0  1 


J 


(353) 


The  vector  31!  is  thus  the  discrete  analog  of  the  second  derivative  of  the 
unknown  vector  Other  forms  for  3  may  be  used,  notably  reference  0*5] . 
As  before,  the  constraint  shall  be  to  minimize  the  quantity  u2(l!)  in  an 
effort  to  achieve  the  smoothest  function  II  with  a  minimum  of  error  in  Zl'.i 

“zM  -  minimum  -g-tf  |  N_  ft  ~  O  (354) 

and  Yz  0  2BTBN  (355) 

::ow,  (340)  and  (355)  in  (354): 

N2  =  (ATA  +  yz5TB)''ATGr  (356) 

This  is  the  required  solution  to  (300)  which  satisfies  the  constraints 
outlined  above.  Again,  setting  ~  6  yields  (305).  A  recent  work 
into  this  particular  solution  has  yielded  considerations  for  the  choice 
of  2^  ,(26]  If  the  factor  26,  is  too  small,  the  instabilities  are  not 
removed  from  the  final  solution.  And  if  chosen  too  large,  the  system 
becomes  overcor strained  and  independent  of  the  measurements.  The  best 


£6 


choice  appears  to  be  /N'  ^ 0-^  and  should  be  so  chosen  until  further 

efforts  provide  a  better  choice. 

Both  solutions  via  tliis  technique  and  the  results  fron  the  next 
method  may  be  improved  in  the  iterative  fashion  as  outlined  in  equations 
(321),  (325)-(327);  IIA.  The  algorithm: 


(initial  estimate) 

1. 

obtain  by  (342),  (356), 

or  (366) 

(m^*1  residual) 

2. 

-  g  -  All^ 

(357) 

(m  error  estimate) 

3. 

^(n)  -  R(n)=>  £(m)_  4-1R(m) 

(353) 

(new  estimate) 

4. 

4(m4l)  _  0(n)  f.(m) 

il  “  T  **# 

(359) 

(test) 


5.  is 


|3WI  <  some  small  positive  number  <f  ? 


no:  m  =  e  -f  1,  go  to  step  2. 

A  (  LJ  ' 

yes:  done,  the  solution  is  I!  '  to  an  accuracy  on  the  order 


of  S. 


The  next  method  involves  the  computation  of  some  statistical  para¬ 
meters.  Gnly  the  general  outline  shall  be  described.  See  reference 
for  details,  let  P,  and  be  positive  definite  matrices  (all  eigen¬ 
values  positive)  that  act  as  weighting  functions,  then: 

vM  =  (G-AN)Tr;CO'/AA/)^C^-No)^  (W-zVo)  (360) 

where  J~,  is  (n)x(n)  and  /Q  is  (k)x(fc). 

Let  IIq  be  the  mean  vector  for  II  as  compiled  from  previous  statistical 
data  (model).  Define  the  expected  value  operator  E  as: 

E  0  s  ^  g6ct)-PC<i)  (361) 

where  f(x)  is  the  probability  density  function  for  g(x).  Let  the  vector 
6  have  zero  mean  and  measure  the  error  in  C-,  Define  the  covariance 
matrices  thus: 

sn=e 


(363) 


=  eM 


Assuming  a  normal  distribution  for  6  and  II,  then  the  Gauss-Karkov 

theorem  yields  a  statistically  optimum  estimate  N  by  minimizing  v(il) 

s 

for: 


then:  s  X  1 AT  “A  tO  * 

where  X  =  +  AT St*  A 

The  resulting  covariance  matrix  for  (II  -  if  )  is  found  to  be  X"1 . 

s 

the  expected  mean  square  error  may  be  written  as: 

E[CM-N,)>-^]=  trace  X'' 


(364) 

(365) 

(366) 

(367) 
Thus, 


(368) 


Without  solving  equation  (300),  one  can  thus  obtain  an  accuracy  estimate 

A 

for  II  under  the  stated  statistical  conditions.  One  can  also  apply  the 

IIA,  (357)  -  (359),  to  enhance  the  estimate  N  , 

s 

The  last  technique  to  be  considered  is  an  iterative  scheme  somewhat 
analogous  to  the  iterative  Improvement  algorithm.  The  method  used  is 
the  Land weber  iteration.  [uj  In  general,  an  iterative  guess,  is 
made  based  on  some  a  priori  knowledge  and  then  improved  successively 
to  form  a  sequence  whose  m^1  member,  U^,  transforms  to  the  measured 
G  as  m  tends  to  infinite.  The  techniques1  success  is  due  to  the  auto¬ 
matic  constraining  of  the  oscillatory  nature  of  for  small  m.  This 


allows  an  appro ximation  to  G  early  enough  to  avoid  the  oscillations  if 
the  sequence  is  terminated  when  the  residual  error  is  on  the  order  of  the 


error  in  G,  The  Landwober  algorithm  uses: 


(369) 


where  III0)  reflects  the  expected  solution.  If  the  coefficient  matrix 


has  all  positive  components,  the  convergence  will  be  improved  by  de- 
* 

fining  a  diagonal  matrix  D: 

ry  .X' 

(370) 


D:  ~  (fr,ai)  j * > 


and  rewriting  the  estimate  as: 

n?1  =  t>flTco-^NVm"1) 


(371) 


This  completes  this  short  survey  of  inversion  techniques.  A  common 
feature  for  all  of  them  is  their  use  of  a  priori  knowledge  to  constrain 
the  selection  of  probable  solutions  in  order  to  enhance  the  choice  of 
the  proper  solution.  Indeed,  this  criterion  must  be  met  by  any  tech¬ 
nique  the  experimenter  chooses  for  inversion,  if  he  is  to  have  any 
confidence  in  the  results. 


V.  conclusion 


Tills  paper  has  endeavored  to  present  the  requirements  for  a  device 
designed  to  map  the  tropospheric  aerosol  properties  on  a  real-time  basis 
The  specific  model  equations  were  derived  after  a  theoretical  devclopmen 
was  performed  in  such  detail  that  the  user  became  familiar  with  the 
limitations  of  the  theory.  A  candidate  for  such  a  device  that  would 
utilise  advantageously  the  specific  theory  characteristics  is  the  laser 
radar,  or  1IDAR.  The  LID Ah  system  allows  for  the  direct  selection  of 

the  key  parameters  of  the  theory  and  yields  as  output  the  required  in- 

« 

puts  for  the  data  inversion  to  begin.  A  few  inversion  techniques  were 
presonted,  but,  duo  to  the  large  available  extent  of  such  methods,  the 
subject  was  by  no  means  exhausted.  The  aerosol  LIDA!  must  be  capable 
of  polarisation  measurements  and  fast  wavelength  scanning  to  increase 
the  number  of  independent  observations  made.  Such  system  capability 
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VS'ST’T^Yr^-  fT‘?'‘=T7?P*^T« — "  'T  T^’r7!  i7  • — n  *r  rr'T?TfnTr>? 

v»wu  u  ^  OrU'.vij-Uw  ii_.U  iLi.'UVili.ilU’i 

This  section  presents  some  further  equations  to  facilitate  the  compu¬ 
tation  of  the  required  fuctions  from  the  Mie  theory. 


Bachs catter:  9-  ISO0 

rrArnl  =  (-/)""  £  <"♦*>  (372) 

T«080°)  =•  J1Ta  Ci8^°)  (373) 

5,  C 190 ° J)  —  n4,  C*0  (n  +  ^)(bn-An)  (371) 

S2(l80’)  =  -5,  (l&0°)  ■  (375) 

cad: scatter  gain:  6  =  @sc.a.O&0  ) -J^rj  j  (-l)  (n*j)(bn-a,^j  (375) 


^7'-) 


1' 


recurrence  relations: 


( nr'.'  \ 


(37?) 


V  *  / 


7n  (cos  e)  r  C35  0  TTn  (C05  0)  -  sm2©^ (cos  e)  (377) 

fr„  Ceos  e)  *  (cos  gX%t)  (“*  e)  -(£l)  K,  Os  e>  (373) 

V^(6ose)=(Zn-i)n-n.,(coia)+ir^tacos  e)  (379) 

1T0  (C05  ©)  =  O  (3C0) 

|T,  (cos  0)  =  !  (321 ) 

t£'(coS0)=O  (3-2) 

'lT;/(C0S©)=  O  (3S3) 

Scattering  coefficients:  The  equations  for  a  and  b  may  ce  written  in 

■  "■  — . . .  n  n  * 

terns  of  the  logarithnic  decrement  function,  (S) .  The  equations 
below  have  been  so  written  and  then  incorporated  into  an  algorithm  that 
was  used  to  compute  the  scattering  ar.d  extinction  efficiencies  for  a 
particular  case  of  c  =  7 •  54(  1— 3.0. C0 1 ).  These  values  were  then  plotted 
on  the  enclosed  graphs  to  picture  the  oscillatory  behaviour  of  the 


|  .Ifyjlfi  ft 


(400) 


’i>* 


.Algorithm.  ,  ,, 

fy(~)  =  (Zn+0  { lanC^^l  +  1**" ^)l  } 

=:>  Q  SCO.  )  ~ 

(J2(<*)  =  ^,(2n*0  [o-nC^)4*  i’nWj 

—>  Qext  W  —  ^2"  ^<?c) 

Begin:  A.  1.  i*?ut:  n  opt  ’  ^ 

2.  compute :  m  x  }  ^  V 
B.  1.  input:  <*><?  (required  accuracy) 

2.  compute:  /?*  ^  &  H 

3.  assign:  n-1  0/  6*0  -  02  ^ 

4.  compute:  s'  6-0  ,  ^  6*) 

x.f")  >  *,c*° 

V?(*) 


1 #  compute 


•h  ft! 


tfC*),  *n  (“> 

3"  00 

.  *«(«) 

.  compute:  £ frt)(V)  *  <J,  («0  +  (2 />*')  +  I  ^  ^  1  ^ 

UPfrO-ftfrOl  ? 

yes:  go  «o  0. l . 
no:  go  to  C*5* 

^1=7 


2.  compute: 

3 

4.  test: 


5.  test:  —S 
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yes:  go  to  D.1 
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